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1. 


INTRODUCTION  . 


The  application  of  the  digital  computer  in  the  field 
of  dance  choreography  is  a  relatively  unexplored  area.  The 
use  of  the  computer  in  other  arts  is  more  widely  known.  The 
computer  has,  for  example,  been  used  to  compose,  arrange  and 
analyze  music,  as  well  as  to  produce  musical  sounds.  It 
has  generated  graphical  works  of  art.  Without  entering  the 
controversy  over  the  computer's  "true  creativity,"  one  can 
state  that  the  computer  has  had  some  impact  on  tho  arts  and 
the  imagination  of  artists.  Its  future  impact  will  no  doubt 
be  greater  as  its  capabilities  become  explored  and  understood. 

In  the  field  of  choreography  the  influence  of  the 
computer  has  so  far  been  negligible.  To  understand  the 
possibility  of  computer  applications  in  this  field,  one 
should  first  consider  the  way  in  which  a  choreographer  works. 

His  modus  ope rand i  is  typically  this:  He  assembles  a 
company  of  dancers,  and  then,  bringing  with  him  only  a  general 
idea  of  what  he  wants,  he  begins  experimenting  and  creating 
his  work  of  art.  The  dancers  are  his  living  tools.  Upon 
completion,  the  record  of  the  evanescent  creation  that  is  a 
dance  lies  only  in  the  minds  of  choreographer  and  dancers 
and,  perhaps,  on  film. 

The  accumulation  of  a  choreographic  literature  has  been 
hampered  primarily  by  the  existence  of  two  interrelated 
problems.  These  are  the  lack  of  an  adequate  movement 
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notation  and  the  enormous  amount  of  information  inherent  in 
the  description  of  even  a  simple  human  movement. 

Dance  notation  systems  do  exist,  the  most  successful 
probably  being  Labanotation ,  devised  by  Rudolf  Laban  and 
first  published  in  1928  [1] .  In  this  system,  positions  of 
various  portions  of  the  body  are  written  on  a  staff  resembling 
a  musical  staff.  The  staff  lines  are  vertical  and  are  read 
from  bottom  to  top.  An  example  of  this  notation  may  be  seen 
in  Appendix  A. 

Transcription  of  dances  into  this  type  of  written 
notation  is  time-consuming  and  difficult,  and  is  not  a  frequent 
practice  of  choreographers. 

Both  of  the  above  mentioned  problems  (the  lack  of  a  good 
notation  system  and  the  large  information  content  of  movement) 
concern  the  manipulation  of  information.  Clearly  here  is 
unbroken  ground  for  the  information  scientist. 

Thus  far,  the  only  significant  use  of  the  computer  that 
has  been  made  in  the  field  of  dance  has  been  to  use  it  to 
produce  lists  of  instructions  for  the  dancer  [2] .  Using  a 
set  of  symbols  representing  possible  movements,  a  computer 
program  generates  a  list  of  movements  to  be  performed  in  the 
given  order.  An  element  of  rdndomness  is  incorporated  into 
the  algorithm  for  producing  these  instructions,  in  an  effort 
to  break  away  from  established  patterns  of  movement.  Such 
instructions  are  quite  general  because  of  the  two  problems 
discussed  above,  thus  leaving  a  large  amount  of  interpretation 
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to  the  dancer. 


Noll  [4]  has  suggested  that  a  different  approach  to 

choreography  be  tried,  making  use  of  the  interactive 

graphical  display  systems  which  are  now  coming  into  use: 

Instead  of  using  the  dancers  as  his  choreographic 
instrument,  the  choreographer  interacts  with  the 
computer  during  the  creative  process.  Stick 
figure  representations  of  the  dancers  appear  in 
some  form  of  three-dimensional  display  on  the  face 
of  an  electronic  display  tube.  The  choreographer, 
by  manipulating  different  buttons  on  the  console 
controls  the  movement  and  progress  of  the  work: 

Different  movements  might  be  stored  in  the  com¬ 
puter's  memory  and  put  together  at  will.  Individual 
movement  restrictions  for  each  dancer  could  even 
be  introduced  into  the  process.  Various  elements 
of  chance  and  randomness  could  be  used  at  the 
discretion  of  ^he  choreographer.  Scenery  and  stage 
props  could  also  be  specified  by  the  choreographer 
and  drawn  by  the  computer.  In  effect,  all  the 
different  aspects  of  the  work  would  be  under  the 
direct  control  of  the  choreographer  to  manipulate 
and  experiment  in  a  myriad  of  combinations.  When 
desired,  certain  portions  of  the  creative  process 
could  be  filled  in  by  the  computer  which  might 
be  programmed  to  learn  the  choreographer's  style. 

All  this  is  a  completely  new  creative  process, 
and  might  result  in  new  dance  forms.  t 

Such  a  tool  would  be  a  solution  to  the  dance  notation 

problem.  It  would  also  be  useful  for  pedagogical  purposes, 

as  a  laboratory  tool  for  student  choreographers,  and  as  a 

device  for  producing  permanent  records  of  dances.  Though 

not  every  choreographer  would  wish  to  use  the  system  described, 

it  would  provide  a  means  of  research  into  human  movement, 

choreographic  styles  and  the  choreographic  process  itself. 

It  is  with  such  a  program  in  mind  that  the  present 

work  was  undertaken.  The  prime  problem  upon  which  attention 

is  focused  is  this:  How  is  a  mathematical  model  of  the 
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human  body  to  be  manipulated  without  resort  to  complex 
mathematical  formulas?  How  can  a  dancer  model  be  controlled 
or  driven  in  a  manner  sufficiently  simple  that  it  not 
interfere  with  the  creative  process  of  the  artist? 

Noll  suggested  manipulating  buttons  on  the  console. 
Another  approach  that  has  been  tried  is  an  anthropometric 
harness  [6).  This  device  is  attached  to  a  person,  and  the 
movements  of  various  key  positions  on  his  body  are  recorded. 
These  movements  are  then  shown  on  an  electronic  display  as 
a  stick  figure  going  through  the  same  motions  as  the  harness- 
wearer.  Possibly  the  interactive  device  known  as  a  joy'stick, 
a  movable  shaft  whose  position  can  be  sampled,  would  be 
useful  in  the  solution  of  the  problem. 

This  dissertation  describes  an  interactive  program 
which  relates  the  angular  movements  of  the  joints  of  the 
display  model  to  hand-drawn  curves.  These  curves  are  drawn 
with  a  stylus  on  a  graphics  tablet.  The  pen  and  tablet  are 
part  of  an  interactive  system  that  includes  a  cathode  ray 
tube  display,  a  digital  computer,  and  the  user  himself.  This 
program  demonstrates  that  this  method  of  model-manipulation 
is  feasible. 

This  dissertation  thus  demonstrates  that  the  notation 
problem  can  be  circumvented  by  the  implementation  of  this 
type  of  graphics  program  as  envisioned  by  Noll.  The  vast 
capacity  of  the  computer  to  store  and  manipulate  large 
amounts  of  complex  information  is  here  brought  to  the  problem 
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of  representing  human  movement.  The  solution  to  this 
problem  is  thus  initiated,  and  it  is  hoped  that  the  effort 
described  here  might  be  of  value  to  future  workers. 
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II.  PROTOTYPE  PROGRAM 


II. 1  Description 

This  program  is  basically  the  type  of  interactive 
program  described  in  the  introduction;  that  is,  a  user  sits 
at  a  display  console  and  views  a  dynamic  display.  The  picture 
in  this  prototype  program  is  a  stick  figure  with  two  moving 
joints.  The  information  that  causes  the  console  to  display 
the  figure  in  its  various  positions  is  generated  in  the 
computer  and  is  based  on  user-supplied  inputs.  These 
inputs  are  primarily  curves  drawn  on  a  tablet  with  a  stylus. 
These  curves  are  shown  on  the  display.  Using  the  coordinates 
of  successive  points  along  the  curves,  the  computer  then 
calculates  arrays  of  angles  through  which  the  moving  parts 

of  the  stick  figure  will  move.  The  illusion  of  movement 

/ 

is  then  created  as  the  computer  displays  in  rapid  succession 
these  various  positions. 

There  are  four  different  basic  displays  that  the  user 
sees.  These  will  be  referred  to  as  "pages"  in  the  remainder 
of  this  discussion. 

The  initial  display  is  page  one,  the  table  of  contents, 
which  merely  lists  the  three  pages  that  follow  it  (Fig.  II. 1-1). 
Pointing  to  the  desired  page  number  or  name  with  the  pen 
causes  that  display  to  appear,  as  if  pages  had  been  turned. 
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Page  two  (Fig.  II. 1-2)  is  concerned  with  determining  the 
movement  of  the  model's  right  leg  about  a  hip  joint.  The 
user  sees  a  small  fixed  figure  with  the  pertinent  joint 
indicated  and  three  rectangles  in  which  the  curves  are  to  be 
drawn.  These  curves  will  represent  the  angular  displacement 
of  the  leg  about  the  X,  Y,  and  Z  axes.  Also  displayed  is  a 
small  diagram  of  these  axes,  showing  their  positive  directions. 
The  remainder  of  the  display  consists  of  the  page  number  and 
a  list  of  user  options.  The  options,  which  are  chosen  by 
pointing  with  the  pen,  include  drawing,  saving  or  erasing 
curves;  returning  to  the  table  of  contents;  and  displaying 
previously  drawn  curves.  Turning  pages  of  the  display  can 
be  accomplished  by  pointing  to  the  arrows  at  the  top  of  the 
picture  by  the  page  number.  The  right  arrow  indicates 
advancing  one  page;  the  left  arrow  indicates  turning  back 
one  page. 
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0*TU  R  NO 

DRAW  CURVES 
KEEP  CURVES 
ERASE  , 

SHOW 

TABLE  OF  CONTENTS 


X  Y 


FIGURE  II. 1-1 
FIRST  DISPLAY 
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« ^rnm+xmt  mrxawm** , 


p  1 
TURN^ 


A  DYNAMIC  BODY  MODEL 

TABLE  OF  CONTENTS 
RIGHT  HIP  ANGLES - 2 

RIGHT  KNEE  ANGLES - 3 

MODEL - 4 

FIGURE  II. 1-2 
SECOND  DISPLAY 
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P  3 

O-TURN-O 

DRAW  CURVES 
KEEP  CURVES 
ERASE 
SHOW 

TABLE  OF  CONTENTS 


MIN  MAX  MIN  MAX  MIN  MAX 


FIGURE  II. 1-3 
THIRD  DISPLAY 
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P  A 
OTURN 


SHOW  MOVE 
ROTATE  X 
ROTATE  Y 
STOP 


TABLE  OF  CONTENTS 


FIGURE  II. 1-4 
FOURTH  DISPLAY 
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DYNAMIC  MODEL  PROGRAM 


INSTRUCTIONS  FOR  THE  USER 


ALL  PAGES: 

*  DEPRESS  PEN  ON  LEFT  OR  RIGHT  ARROW  TO  TURN  PAGE  BACK  OR  FORWARD. 

*  DEPRESS  PEN  ON  "TABLE  OF  CONTENTS"  TO  RETURN  TO  PAGE  ONE. 

PAGE  1 

*  DEPRESS  PEN  ON  NAME  OR  NUMBER  OF  PAGE  DESIRED. 

PAGE  2 

DEPRESS  PEN  ON  DESIRED  OPTION: 

- ►*  "DRAW  CURVES";  THEN  DRAW  CURVES  IN  BOXES  IN  ANY  ORDER.  IT 

IS  NOT  NECESSARY  TO  ERASE  BEFORE  CHOOSING  THIS  OPTION.  ONE 
OR  MORE  CURVES  MAY  BE  RE-DRAWN;  THEN  CHOOSE  "KEEP  CURVES." 

*  "KEEP  CURVES"  IF  CURVES  DRAWN  ARE  ACCEPTABLE. 

*  "ERASE"  REMOVES  ANY  CURVES  DRAWN. 

*  "SHOW"  DISPLAYS  LAST  CURVES  DRAWN.  THIS  IS  USED  WHEN 
RETURNING  FROM  ANOTHER  PAGE  OR  AFTER  ERASURE. 

*  "RE-USE"  THE  CURVES  KEPT  BEFORE  WILL  BE  USED  WITH  DIFFERENT 
CURVES  DRAWN  ON  ANOTHER  PAGE. 

PAGE  3 

DEPRESS  PEN  ON  DESIRED  OPTION: 

*  "DRAW  CURVES";  THEN  DRAW  CURVES  IN  THE  BOXES,  ANY  ORDER. 

*  "KEEP  CURVES"  IF  CURVES  JUST  DRAWN  ARE  ACCEPTABLE. 

*  "ERASE"  REMOVES  ANY  CURVES  DRAWN  FROM  THE  DISPLAY. 

*  "SHOW"  DISPLAYS  LAST  CURVES  DRAWN. 

*  "RE-USE"  THE  CURVES  KEPT  BEFORE  WILL  BE  REPROCESSED  AFTER 
SOME  CHANGE  MADE  ELSEWHERE. 


PAGE  4 

DEPRESS  PEN  ON  DESIRED  OPTION: 

*  "SHOW  MOVE"  DISPLAYS  MOVING  FIGURE. 

*_*  "ROTATE  X"  MOVES  VIEWPOINT  UP  OR  DOWN  ACCORDING  TO  VALUE  TYPED. « - 

*-*  "ROTATE  Y"  MOVES  VIEWPOINT  LEFT  OR  RIGHT  ACCORDING  TO  VALUE  TYPED.* 

*  "STOP"  TERMINATES  PROGRAM. 

"RE-USE"  COMPUTES  NEW  VIEW  AFTER  SOME  CHANGE  BELOW. 

TELETYPE  OPTIONS;  AFTER  TYPING  HIT  RETURN  KEY 
I-*  EYESTG-n  CHANGES  DISTANCE  EYE  TO  STAGE,  WHERE  n  IS  IN  FEET. 

INITIAL  SETTING  IS  20  FEET. 

V-*  EYETUB=n  CHANGES  DISTANCE  EYE  TO  DISPLAY  TUBE.  INITIALLY 
1.4  SCREEN  DIAMETERS. 

*  LEFT=n,  WHERE  n  IS  BETWEEN  -.1  and  .42 - 1 - •* 

MOVES  EYE  RIGHT  OR  LEFT.  INITIAL  SETTING  0. 

*  UP=n,  WHERE  n  IS  BETWEEN  -.32  and  .02 - - - 

MOVES  EYE  DOWN  OR  UP.  INITIAL  SETTING  0. 

► — *  NN*=n  CHANGES  SPEED  OF  MOVEMENT,  WHERE  n  IS  BETWEEN  4  and  25, 

THE  LARGER  NUMBER  BEING  SLOWER.  INITIALLY  25. 


USER'S  MANUAL 
FIGURE  II. 1-5 
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Page  three  (Fig.  II. 1-3)  is  similar  to  page  two,  but 
is  concerned  with  the  angle  at  the  knee.  As  movable  joints 
are  added  to  the  program,  similar  pages  may  be  added.  The 
program  is  so  structured  that  it  may  be  extended  in  this  way. 
This  will  be  discussed  further  in  Section  III. 

Page  four  (Fig.  II. 1-4)  displays  a  larger  stick  figure 
and  a  set  of  user  options .  If  the  option  to  display  the 
current  movement  is  selected,  the  stick  figure  is  shown  in 
varying  positions  in  rapid  succession,  so  that  he  appears  to 
move.  This  dynamic  display  may,  of  course,  be  repeated  as 
many  times  as  desired.  Page  four  has  the  usual  options  for 
changing  the  display  to  other  pages. 

Also  associated  with  page  four  are  options  to  change 
the  positions  of  the  viewpoint  and  the  stage,  and  to  alter 
the  speed  of  the  figure's  movement.  These  will  be  explained 
more  fully  in  Sections  II. 2-4  and  II. 4. 

As  he  sits  at  the  console#  the  user  has  beside  him  a 
one-page  manual  of  instructions  for  using  the  program.  This 
manual  is  intended  for  use  by  non-technical  persons  with  no 
programming  experience,  and  is  shown  in  Figure  II. 1-3. 

II. 2  Data  Definition  and  Manipulation 

All  objects  in  this  program  are  ultimately  described 
in  a  three-dimensional  coordinate  system  having  its  origin 
at  the  viewpoint.  The  stage  and  each  joint  of  the  stick 
figure  have  associated  local  coordinate  systems.  For  the 
hip  joint,  this  local  system  originates  at  the  base  of  the 
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spine,  and  the  vertical  (Y)  axis  corresponds  to  the  spine. 
For  the  knee  joint,  the  system  has  its  origin  at  the  knee 
point,  and  the  vertical  axis  corresponds  to  the  upper  leg 
vector ,  whatever  orientatioil  that  vector  may  have. 

In  addition,  each  joint  has  associated  with  it  the 
.following : 

1  initial  point  coordinate  vector 
3  x  n  array  of  angles 
6  angle  constraints 

3x3  natural  linkage  rotation  matrix 
4x4  update  transformation  matrix 

3  table  pointers 

4  x  n  display  point  array 

where  n  is  the  number  of  points  in  time  for  which  sets  of 
display  points  have  been  computed.  These  will  be  explained 

I 

in  the  sections  that  follow. 

II. 2-1  Angle  Determination 

Each  moving  joint  has  associated  with  it  an  array  of 
angle  values,  representing  the  different  positions  of  the 
limb  or  vector  at  different  points  in  time.  These  angle 
values  are  relative  to  the  local  coordinate  system.  Each 
angle  is  divided  into  three  components  representing  the 
angular  displacements  about  the  three  axes.  Each  of  these 
three  components  may  be  thought  of  as  an  array  corresponding 
to  one  of  the  user-drawn  input  curves. 

The  curve-drawing  rectangles  ara  labeled  "X,"  "Y,"  and 
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I 

"Z"  for  the  axes  about  which  the  angular  displacements  are 
measured.  The  vertical  arrows  seen  on  the  display  indicate 

t 

that  the  positive  direction  of  the  independent  variable, 
time,  is  from  the  bottom  to  the  top  of  the  rectangle.  This 
orientation  was  chosen  to  correspond  to  that  of  Labanotation , 
which  is  also  read  from  the  bottom  upwards.  Within  each 
rectangle  is  a  vertical  dashed  line,  indicating  the  reference 
angle  of  zero  degrees  of  rotation  about  that  axis.  The  left 
and  right  sides  of  the  rectangle,  labeled  minimum  and  maximum, 
represent  angles  of  100%  of  the  constraint  angles. 

A  constraint  angle  is  set  by  human  limitations  and  is 
the  limiting  angular  displacement  of  a  limb  about  an  axis. 

In  this  program  each  angular  component  has  negative  and 
positive  constraints.  The  left  (user's  left)  side  of  the 
displayed  rectangle  represents  the  constraint  when  the  limb 
is  rotated  in  a  negative  direction  about  that  axis;  the  right 
side,  positive.  Some  joints,  such  as  the  knee,  will  have  at 
least  one  constraint  of  zero  degrees. 

Thus  it  is  the  curve's  amplitude,  or  horizontal  distance 
of  the  curve  from  the  dashed  line,  that  determines  the  angular 
component  at  each  time  point.  This  component  represents 
the  total  angular  displacement  about  one  axis  at  one  time 
point  and  is  made  to  be  directly  proportional  to  the  input 
curve  amplitude. 

The  angular  component  is  computed  from  this  formula: 

0  =  (a/d) K 

where  0  is  the  angular  component,  a  is  the  amplitude  of  the 
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input  curve  (in  display  scope  units) ,  d  is  the  distance 
between  the  dotted  line  and  box  edge  (in  scope  units) ,  and 
K  is  the  constraint  angle' (in  radians)  in  either  positive  or 
negative  rotation,  depending  on  which  side  of  the  dashed  line 
the  curve  point  is. 

The  angle  constraints  are  constants  entered  in  just  one 
place  in  the  program,  in  the  data  section,  and  are  easily 
altered. 

The  dashed  line  is  located  to  reflect  accurately  the 
ratio  between  the  absolute  values  of  negative  and  positive 
constraints.  If  constraints  are  changed,  the  dashed  line  is 
automatically  repositioned. 

It  is  seen  that  this  method  of  angle  determination 
obviates  the  need  .for  checking  angles  computed  against 
limiting  constraints.  It  is  unnecessary  to  guard  against 
computing  unrealistic  positions,  because  it  is  impossible 
to  exceed  the  constraints.  Constraints  will  be  discussed 
further  in  section  III. 

II. 2-2  Vector  Determination 

Homogeneous  coordinate  representation  is  used  to  repre¬ 
sent  point  coordinate  vectors  and  transformation  matrices 
[7]  [8].  Homogeneous  representation  is  simply  the  repre¬ 
sentation  of  an  n-dimensional  object  in  n+1  coordinates.  For 
example,  a  point  coordinate  in  three-dimensional  space  is 
represented  by  four  coordinates,  the  fourth  one  being  the 
homogeneous  coordinate . 
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The  reason  for  using  homogeneous  coordinates  is  that 
all  rotations  and  translations  of  a  particular  vector  can 
be  represented  in  a  single  matrix. 

Each  limb  is  represented  initially  as  a  row  vector  in 
local  coordinates,  and  is  oriented  correctly  by  multiplying 
it  by  its  updated  transformation  matrix,  which  contains 
all  the  rotations  and  translations  through  which  it  has 

passed,  including  that  of  previous  systems  to  which  it  is 

) 

related. 

The  transformation  matrix  can  be  viewed  as  consisting  of 
components : 

|3  x  3 

I  Rotation 
Matrix 

X  Y  Z 


0 

0 

0 

1 


The  upper  left  3x3  matrix  represents  the  rotation,  while  the 
first  three  values  in  the  bottom  row  represent  the  trans¬ 
lation  in  the  X,  Y,  and  Z  directions,  respectively.  (The 
last  column  represents  perspective  transformations  but  is  not 
used  as  such  in  this  work  [10]  .) 

The  rotation  matrices  used  were 


about  the  X  axis : 


about  the  Y  axis: 


'  1  0  0* 

0  C  S 

.0  -S  C 

'  C  0  -S 

0  10 

6  0  C 
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about  the  Z  axis:  C  S  0 

-SCO 
.  0  0  1 

where  C  =  cos  0,  S  =  sin  0.  0  =  angular  displacement  about 

one  axis. 

An  example  of  a  unit  vector  on  the  X  axis  being  rotated 
about  the  Y  axis  follows. 

U  0  0  1]  C  S  0  0  =  [C  SOI] 

-S  C  0  0 

0  0  10 

0  0  0  1. 

Translating  this  resultant  vector  a  units  along  the  X  axis, 
b  along  the  Y  and  d  along  Z  would  be  performed  as  follows. 

[C  SOI]  1000  =  [C+a , S+b ,d , 1] 

0  10  0 

0  0  10 

.a  d  1  j 

Both  operations  could  be  performed  by  one  multiplication, 
where  the  transformation  matrix,  M,  is  the  product  of  R,  the 
4x4  rotation  matrix  (in  homogeneous  representation) ,  and  T, 
the  translation  matrix: 

M  :=  R  x  T 
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For  the  above  example , 


C  S  0 


M 


C 

0 


0 

1 


a  b  d 


0 

0 

0 

1 


and 


[10  0  1] 


C  S  0  0 

-S  C  0  0 


=  [C+a , S+b ,  d ,  1  ] 


0  0  10 

a  b  d  1 


as  before. 

Note  that  this  represents  the  rotation  first  followed  by 
the  translation.  The  reverse  order  would  give  a  different 
result. 

The  rotation  matrix  is  the  product  of  the  natural  linkage 
rotation  matrix,  the  three  rotational  components  of  the  local 
system,  and  the  rotation  matrix  of  the  previous  or  "father" 
system  to  which  the  local  system  is  linked.  For  example,  if 
the  local  system  is  that  one  with  its  origin  at  the  knee,  the 
rotation  matrix  of  the  hip  joint  must  be  included  in  the 
product  that  forms  the  knee  updated  rotation  matrix.  In  turn, 
the  rotation  matrix  of  the  hip  joint  includes  rotations  of  the 
spine  relative  to  the  stage,  and  of  the  stage  relative  to 
the  viewpoint. 
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The  natural  linkage  rotation  matrix  describes  the 
orientation  of  the  local  coordinate  system  initially  (at 
rest)  to  the  system  just  preceding  it.  These  are  all 
identity  matrices  in  this  program,  but  they  are  included  for 
generality. 

The  translational  portion  of  the  transformation  matrix 
represents  the  origin  of  the  local  coordinate  system. 

Once  the  transformation  matrix  is  set  up,  the  position 
of  the  vector  it  represents  is  determined  by  the  following 
multiplication , 

Vk  =  VD  X  [M] 

where  V  is  the  initial  point  coordinate  row  vector,  M  is  a 
4x4  transformation  matrix,  and  is  the  translated  and 
rotated  point  coordinate  vector.  Vq  is  represented  in  local 
coordinates;  for  example, 

Vq  =  [100,  0,  -5,  1] 

would  represent  a  point  that  was  translated  100  scope  units 
from  the  origin  along  the  X  axis,  zero  along  the  Y  axis,  and 
minus  five  along  the  Z  axis  of  its  associated  local  coordinate 
system.  The  homogeneous  coordinate  in  its  most  general  case 
is  a  scale  factor  by  which  the  entire  vector  is  divided,  but 
in  this  work  it  will  always  have  a  value  of  one. 

11.2-3  Perspective  Transformation 

The  points  obtained  by  the  process  just  described  have 
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coordinates  based  on. a  system  with  its  origin  at  the  eye  or 
viewpoint. .  These  are  converted  to  display  points  for  a 
perspective  picture  [11]  with  the  following  formulas: 


-  xs 

Y-a  -  (Yv  Zs)/Zj  -_YS 


where  ,  Y^)  is  the. -displayed  point  in  screen  coordinates, 

Yjj,  Zj.),  is  the  object  point  in  eye  coordinates,  and 

(X  ,,  Y.  ,,  Z:  )  is  the  screen  origin  in  eye  coordinates  (See 
s  s.  s_ 

Eig..  IL.  2.,  .3  -1}  . . 

The.  values,  of.  Z  and  Z  .  (the  distances  from  viewpoint 

Sf  3 

to  screen. and  object)  may  be  varied  while  using  the  program, 
as- will  be  explained  in  the  next  section.  The  user  may  thus 
obtain,  the  perspective  from  the  distance  desired,  and  experi¬ 
ment  wifcn  subjective  effects  obtained  by  varying  the  viewing 
distance.. 


II., 2- 4;  Camera.  Positioning 

The  ultimate  origin  of  all  points  in  all  coordinate 
systems  is  the  viewpoint.,  The  viewpoint,  therefore,  cannot 
be  altered..  The  stage,,  however ,  may  be  repositioned  relative 
to  th.e.  viewpoint,,  which,  creates  the  illusion  of  moving  the 
camera  or  viewpoint.. 

There  are  four  parameters  which  may  be  changed  to  vary 
the  positions  of  the  stage  and  of  the  display  scope  relative 
to  the  viewpoint..  These  must  be  typed  in  on  the  Teletype 
while  page  four  is  being  displayed.  Precise  instructions  for 
doing  this  are  to  be  found  in  the  user's  manual. 

\\ 


The  first  of  these  parameters  is  the  distance  from  eye 
to  stage.  It  is  represented  in  feet,  and  is  initially  set 
to  twenty  feet. 

The  distance  from  eye  to  display  tube  may  likewise  be 
altered  by  typing  in  a  new  value.  This  value  is  initially 
set  at  1.4  x  width  of  display  tube. 

The  stage  may  be  rotated  about  the  X  or  Y  axes  of  the 
eye's  coordinate  system.  These  rotation  angles  are  in 
radians  and  are  initially  set  to  zero.  Suggested  limits 
within  which  to  contain  the  rotation  angles  are  included  in 
the  user's  manual.  Since  they  are  represented  in  the  same 
form  as  limb  rotation  angles,  they  are  actually  angle  arrays. 
The  data  structure  thus  allows  for  panorama  viewing,  though 
this  program  implements  only  a  static  viewpoint. 

II. 3  Data  Structure 

Each  of  the  origins  of  the  various  coordinate  systems 
is  associated  with  a  node.  The  nodes  are  related  to  each 
other  in  a  tree  structure  (Fig.  II. 3-1)  and  are  triply 
linked  by  father-brother-son  pointers  as  shown  in  Table 
***3-l.  Each  node  except  the  first  points  to  a  father  node, 
to  one  "younger"  brother  and  to  a  first  son.  The  pointer 
value  of  zero  in  the  father  column  indicates  that  that  node 
is  the  origin  of  the  tree.  Pointer  values  of  "  -  "  in 
brother  and  con  columns  indicate  that  that  relationship  is 
unfilled ,  and  they  are  actually  represented  in  the  program 
by  values  of  zero. 
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ALGORITHM  FOR  DETERMINING  COORDINATES  OF  DISPLAY  POINTS 
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The  algorithms  for  determining  the  display  points  and 
for  displaying  them  are  both  driven  by  this  table.  The 
point-determinatior.  algorithm  is  shown  in  Figure  II.  3-2, 
and  the  display  algorithm  is  similar.  A  more  complex  display 
could  be  handled  by  these  key  routines  by  entering  the  appro¬ 
priate  table  and  initial  data  in  the  data  section  of  the 
program. 

II. 4  Curve-Processing 

When  a  curve  is  drawn  on  the  tablet  with  the  stylus, 
the  information  is  stored  as  an  array  of  three-dimensional 
(X,  Y,  Z)  points.  The  value  of  the  Z  coordinate  is  infor¬ 
mation  about  the  continuity  {beginning  and  ending)  of  the 
curve  and  is  discarded  after  use.  The  X  and  Y  coordinates 
are  the  usual  horizontal  and  vertical  coordinates  (in  scope 
units) . 

Once  a  point  is  recognized,  a  succeeding  point  is  not 
recognized  until  the  stylus  has  moved  outside  of  a  4  x  4 
square  area  centered  at  the  preceding  point. 

It  will  be  seen,  then,  that  the  number  of  points 
recognized  on  one  of  the  user-drawn  curves  will  vary  from 
curve  to  curve.  A  straight  line  will  consist  of  fewer 
points  than  a  curve  with  many  inflections.  The  problem 
thus  arises  of  relating  the  curve  points  to  each  other. 

It  is  desired  that  all  curves  contain  the  same  number 
of  points  and,  further,  that  points  with  like  vertical  coordi¬ 
nates  will  correspond  to  each  other.  This  is  accomplished 

26 


by  thinning  the  points  to  a  constant  number  of  points  which 
are  evenly  spaced  along  the  vertical  axis.  Points  are  then 
filled  in  between  these  thinned  points  in  the  following  . 
manner:  A  parabolic  approximation  tp  the  first  four  points 
is  obtained  by  the  least-squares  method  [12],  [13].  Five 

points  are  then  interpolated  in. the  middle  interval,  and 

\ 

the  procedure  is  repeated,  for  all  succeeding  intervals,  save 
the  last.  Interpolation  is  linear  inn the  end  intervals. 

The  number  of  points  to  which  the  curve  is  thinned  is 
set  to  twenty-five  in  the  program, ;  but- can  be  changed  by 
typing  in  different  values  on.  the  Teletype.  (Instructions 
for  so  doing  are  in  the  user's  manual-).  Changing  this 
parameter  alters  the  number  of  "frames^'”  shown  when  the 
movement  is  displayed.  This  parameter.- i Sj  then,  in  effecp 
a  scaling  factor  which  varies  the  spe^d  of.  the  motion. 


11 


III.  AUGMENTING  THE  PROTOTYPE 


The  prototype  program  is  intended  as  a  suggested  approach 
to  the  solution  of  the  problem  of  manipulating  the  figures  in 
a  choreography  program.  There  are  many  factors  to  be  con¬ 
sidered  in  further  development  of  this  idea.  Some  of  these 
factors  will  be  discussed  in  this  section. 

III.l  Improving  the  Existing  Model 

There  are  a  number  of  significant  improvements  that 
could  be  made  within  the  context  of  the  existing  12-node  data 
structure  and  the  stick  figure  representation  of  the  model 
(See  Figure  II. 3-1). 

Full  articulation  of  the  existing  model  is  desirable, 
and  the  existing  structure  is  such  that  achieving  this  is 
primarily  a  clerical  effort  on  the  part  of  the  programmer. 

Two  things  must  be  added  to  the  program  for  each  additional 
joint  articulated:  angle  constraints  and  a  new  display 
(curve-drawing  page) ,  analogous  in  every  way  to  the  second 
display  (page  two) .  The  angle  constraints  are  entered  in 
the  Block  Data  section  of  the  program.  All  necessary  arrays 
are  already  in  existence  and  correctly  dimensioned  for  full 
articulation. 

Implicit  in  full  articulation  is  the  movement  of  the 
1  # 

dancer  relative  to  the  stage.  It  will  be  recalled  that 
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there  is  a  separate  stage  coordinate  system  which  is  between 
the  eye  (viewpoint)  coordinate  system  and  the  initial  or 
base-of-spine  coordinate  system  of  the  dancer.  Therefore, 
movement  of  the  dancer  across  the  stage  is  possible  within 
the  data  structure,  though  it  is  not  implemented  in  the 
prototype . 

Such  movement  would  be  internally  represented  in  trans¬ 
formation  matrices,  as  are  all  other  movements.  In  addition, 
a  special  array  of  vectors  would  be  needed  to  represent  the 
initial  point  coordinate  vector,  since  this  is  not  of  constant 
magnitude.  The  problem  arises  of  entering  the  latter  type 
of  information  in  the  program.  How  does  the  user  indicate 
the  desired  movement  of  the  whole  dancer  relative  to  the 
stage? 

Copeland  [14]  has  written  a  graphics  program  which 
traces  the  path  of  the  stylus  on  a  tablet  area  representing 
the  stage  in  plan  view;  one  or  more  small  triangles,  repre¬ 
senting  dancers,  then  traverse  this  path.  Such  a  program 
could  be  used  to  put  in  information  about  the  X  and  Z 
translation  (upstage/downstage,  right/left),  and  the  Y 
translation  (vertical  movement)  could  be  found  in  a  similar 
manner.  A  hand-drawn  curve  on  the  tablet  could  be  interpreted 
as  a  simple  linear  relationship  between  time  and  the  distance 
of  the  base  of  the  spine  from  the  floor.  Obviously,  rather 
severe  constraints  would  have  to  be  built  in  here.  Perhaps 
the  program  could  check  the  feasibility  of  these  values  for 
vertical  movement  by  comparing  them  with  the  corresponding 
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leg  positions.  Eventually,  such  factors  as  momentum, 
gravity,  and  muscular  strength  might  be  included  in  this 
feasibility  check. 

Storage  requirements  of  an  expanded  program  will  need  to 
be  watched.  The  prototype  program  requires  a  little  less 

t 

than  2,1,000  words  of  storage.  This  includes  the  expanded 
dimensions  mentioned  before.  The  Univac  1108  has  a  core 
storage  capacity  of  65,000  words  of  which  approximately 
53,000  are  available  to  the  user.  Each  new  display  (for 
each  new  joint  articulated)  adds  about  3,700  words  to  the 
storage  requirement,  about  half  of  which  is  the  new  page 
display  for  curve-drawing,  and  half  of  which  is  code.  Thus, 
any  eight  of  the  existing  unarticulated  nodes  could  be  articu¬ 
lated  without  requiring  additional  storage,  but  articulation 
of  all  twelve  nodes  would  not  be  possible  without  some  space¬ 
saving  changes. 

Some  storage  economy  could  be  effected  by  packing  the 
object  point  coordinates  two  or  three  to  a  (36-bit)  word. 

This  is  possible,  because  these  coordinates  are  represented 
internally  as  signed  integers  with  known  ranges.  A  limit 
would  need  to  be  placed  on  the  Z-coordirjate ,  the  distance  from 
viewpoint  to  figure.  (It  will  be  recalled  that  the  user  may 
vary  this  distance.)  If  it  were  limited  to  fifty-eight  feet, 
the  coordinates  could  be  packed  three  to  a  word.  Packing  two 
to  a  word  would  permit  a  limit  of  3745  feet. 

Another  improvement  that  could  be  made  in  the  existing 
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program  is  to  make  the  input  language  entirely  stylus-and- 
tablet  oriented.  The  need  for  this  becomes  obvious  to  a  user 
attempting  to  exercise  the  Teletype-input  options,  such  as  by 
varying  the  distance  of  the  stage.  Though  the  Teletype  is 
located  immediately  next  to  the  tablet  and  the  user,  it  is 
disconcerting  to  alternate  between  input  devices.  Programming 
this  change  would  not  be  difficult. 

Another  change  concerns  matrix  multiplication.  When  a 
rotation  matrix  undergoes  a  number  of  computer  multiplications, 
as  it  doos  in  this  program,  the  orthogonality  properties  are 
lost  because  of  computer  truncation.  Mahl  [15]  has  devised 
a  correction  equation  which  restores  orthogonality  to  t.he 
rotation  matrix.  It  is  a  rapidly  converging,  iterative  scheme 
based  on  a  Taylor's  series.  As  the  numoer  of  matrix  multi¬ 
plications  is  increased  in  expanding  the  model,  incorporation 
of  this  scheme  into  the  program  might  be  considered. 

III. 2  Anthropometric  Considerations 

A  truer  representation  of  the  mechanical  relationships 
of  the  human  body,  or  in  other  words,  a  better  stick  figure, 
should  be  an  early  goal  of  one  improving  the  prototype. 

III. 2-1  Constraints  on  Movement  of  Joints 

In  the  prototype  program,  crude  approximations  of  the 
natural  constraints  on  joint  movements  are  used,  the  model 
being  a  crude  approximation  of  a  human  figure. 

These  limits  on  joint  movements  were  assumed  to  be 
constant,  although  in  reality  most  constraints  will  vary 
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according  to  the  position  of  other  parts  of  the  body.  The 
constraints ' will  of  course  also  vary  from  individual  to 
individual,  and  even  from  day  to  day  within  the  same  individual 
[17] .  However,  the  working  assumption  of  constant  constraints 
is  probably  no  grave  drawback  even  in  a  highly  developed 
program  if  the  constraints  supplied  are  quite  liberal.  The 
end  product  of  the  program  is  a  person's  subjective  impression 
of  a  graphical  display.  Unsatisfactory  postures  can,  hope¬ 
fully,  be  altered  by  the  user,  until  they  appear  to  be  satis¬ 
factory  from  a  variety  of  vantage  points . 

7 

These  angular  constraints  for  joint  movement  have  been 
extensively  measured  for  medical  purposes  and  for  use  in 
industrial  and  prosthetic  design.  Examples  are  to  be  found 
in  [16],  [17],  [18],  and  [19].  The  last  two  give  average 
normal  ranges,  which  an  average  individual  could  be  expected 
to  satisfy.  A  dancer  would  have  less  restrictive  constraints. 
Reference  [16]  gives  the  2.5  percentile,  50  percentile  and 
97.5  percentile  ranges  for  men,  for  women,  and  for  children 
of  four  different  ages.  This  reference  also  has  a  good 
bibliography  of  anthropometric  studies.  Reference  [17]  also 
gives  constraint  data  for  different  percentiles  of  the  popu¬ 
lation  studied. 

III. 2-2  Natural  Linkages 

It  will  be  recalled  that  the  prototype  model  has  a 
coordinate  system  associated  with  each  joint.  When  the 
figure  is  at  rest,  the  various  rotation  angles  are  zero. 
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However,  this  resting  position,  or  position  of  no  rotation, 
may  orient  the  local  coordinate  system  at  an  angle  with 
respect  to  the  system  with  which  it  is  linked.  The  natural 
orientation  or  at-rest  rotation,  which  is  assumed  to  be  zero 
degrees  in  the  prototype,  will  need  to  be  determined  and 
entered  as  data  in  an  improved  program.  This  writer  was 
unable  to  find  such  data  in  the  literature. 

When  the  dancer  representation  has  reached  a  level  of 
sophistication  such  that  the  above-mentioned  data  is  being 
incorporated  into  the  program,  it  would  be  appropriate  to 
experiment  with  permutations  of  the  order  of  multiplication 
of  the  three  matrices  that  make  up  any  individual  rotation 
matrix.  This  order  remains  fixed  in  the  prototype  as  the 
rotation  about  the  X  axis  first,  followed  by  the  rotations 
about  Y  and  the  Z.  However,  it  might  be  discovered,  for  k» 
example,  that  the  movement  of  the  wrist  is  most  naturally 
described  by  rotating  first  about  the  Z  axis,  followed  by 
X  and  then  Y.  (There  are  six  possible  permutations  of 
the  three  rotations.)  Note  that,  matrix  multiplication 
not  being  commutative,  the  product  or  result  will  vary, 
depending  on  the  order  in  which  the  rotations  take  place. 

III. 2-3  Body  Dimensions 

Data  for  such  things  as  limb  length,  shoulder  width, 
etc.,  may  also  be  found  in  anthropometric  literature  [16], 
[17] .  (Limb  length  is  not  truly  fixed  as  an  arm  or  leg 
moves,  but  again,  it  may  be  assumed  constant,  probably 
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without  much  loss  of  verisimilitude,  particularly  in  the 
context  of  currently  available  display  resolution.) 

Using  these  kinds  of  measurements,  a  mathematical 
representation  of  an  improved  stick  figure  could  be  created. 

Or,  going  one  step  further,  a  "f leshed-out "  figure  could  be 
represented,  still  using  only  straight  lines. 

It  is  possible  within  the  general  framework  of  the 
program  to  create  a  wire-frame  or  a  polyhedral  representation 
of  the  dancer.  Either  is  envisioned  as  a  set  of  rigid,  non¬ 
intersecting  three-dimensional  structures.  The  single  vectors 
of  the  prototype,  each  of  which  is  associated  with  a  different 
coordinate  system,  would  each  in  the  expanded  version  be  a  set 
of  vectors,  all  fixed  in  one  coordinate  system.  As  soon  as 
the  algorithm  for  determining  point  coordinates  has.  found  the 
primary  point  coordinate  vector,  all  the  other  point  coordinate 
vectors  fixed  in  that  same  coordinate  system  should  be 
determined.  (See  Figure  II. 3-2.)  The  display  algorithm  would 
likewise  have  to  be  altered  by  including  for  each  coordinate 
system  a  table  of  the  mapping  of  the  points  associated  with 
that  system,  which  mapping  is,  in  the  current  program,  a 
single  vector.  A  wire-frame  or  a  polyhedral  figure  could  be 
described,  depending  on  the  mapping  used.  Included  in  this 
mapping  could  be  inter-limb  connections  to  create  knees, 
elbows t  etc.  The  vectors  or  surfaces  that  comprise  these 
connections  are,  in  effect,  flexible  coverings  of  skin. 

The  curtains  and  scenery  (stage  props)  description 
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could  be  included  within  this  expanded  data  structure.  The 
various  elements  of  stage  and  scenery  display  would  be 
vectors  fixed  in  the  stage  coordinate  system,  corresponding 
co  the  set  of  vectors  belonging  to  a  single  limb. 

111. 3  Display  Considerations 

Along  with  wire-frame  representation  of  the  dancer, 
hidden-line  removal  would  be  desirable,  so  that  the  dancing 
figure  would  be  displayed  not  as  a  transparent  wire-frame 
structure  but  rather  as  a  solid  structure.  Alternatively, 
if  the  object  points  were  mapped  in  a  polyhedral  represen¬ 
tation,  hidden  surface  removal  would  be  appropriate.  Either 
of  these  approaches,  as  well  as  shaded  renderings,  is  possible 
with  existing  algorithms. 

Another  feature  that  could  be  added  to  the  prototype 
program  is  a  provision  for  windowing.  Because  of  the  wrap¬ 
around  characteristic  of  the  available  display  hardware, 
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the  model  will,  when  moved  off  the  display,  reappear  on  the 
opposite  side.  Thus  there  is  a  need  to  use  some  existing 
hardware  or  software  which  will  perform  windowing  or 
clipping;  that  is,  tha-  will  eliminate  portions  of  the  picture 
outside  of  the  viewing  window  [11] . 

111. 4  Library-based  Program 

Ultimately,  the  type  of  program  under  consideration 
would  have  to  be  based  on  a  library  of  standard  movements 
if  it  is  to  be  a  useful  program.  The  user  would  select  the 
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desired  standard  moves,  perhaps  modify  them  in  some  way, 
and  then  join  them.  The  user  would  also  have  the  option  of 
creating  new  moves,  of  course,  perhaps  adding  them  to  his 
library.  Maximum  flexibility  would  have  to  be  allowed  by 
providing  options  for  timing,  splicing,  repeating,  saving, 
altering  and  interpolation  between  movements  by  the 
p  rogram. 

There  should  be  provided  a  choice  of  different  dancer 
models,  perhaps  allowing  the  choreographer  to  specify 
individual  body  parameters.  Different  standard  scenery 

elements  or  the  option  of  uniquely  specifying  them  should 

( 

likewise  be  available. 

The  standard  library  movements  could  be  created  in  a 
number  of  ways,  and  it  remains  for  future  workers  to  determine 
the  best  approach.  Possibilities  include  the  curve-drawing 
approach  of  this  paper,  the  anthropometric  harness  [6],  and 
the  joystick.  The  track  ball  is  another  input  device  which 
might  be  considered,  although  the  author  found  it  unsatis¬ 
factory  in  a  preliminary  trial.  Perhaps  analysis  of  multiple 
photographs  of  live,  moving  dancers  will  provide  the  best 
approach  to  compiling  a  library  of  movements. 

The  method  or  methods  of  joining  the  standard  movements 
might  be  different  from  the  methods  used  in  creating  them. 

The  possibilities  include,  in  addition  to  those  mentioned  in 
the  previous  paragraph,  automatic  filling-in  and  joining  by- 
the  program  and  mathematical  specifications  of  the  user. 

This  problem  of  joining  or  splicing  was  not  studied  per  se 
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in:,  this  work,  but  the  curve-drawing  approach  provides  a  tool 
for  accomplishing  this.. 

The  implementation  of  a. graphics^ program  that  would  be 
truly-  useful  to  choreographers  is  a  profound  undertaking  that 
has:  not  yet-  been • achieved . .  However  , ,  computer  applications, 
have,  been  thus  far  limited  more  by  the^ imagination. of  men-, 
than'. by/ the  computer:- hardware ,  and  the-- realization  of  a 
useful  choreography,  program  of  the  type  outlined  may  be 
realistically - anticipated . . 
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APPENDIX  A 


An  example  of  Labanotation  is  shown.  This  is  taken  from 
page  143  of  reference  [1] .  It  describes  one  dancer's 
movements  for  several  seconds. 

This  is  reprinted  with  the  permission  of  the  pub¬ 
lishers  ,  Theatre  Arts  Books,  New  York,  copyright  1954  by 
The. Dance  Notation  Bureau,  Inc. 

STUDY  IN  BALLET  SI  YLE 
’  (Use  of  Parts  of  the  Leg) 

Suggested  music;  The  Coda  of  Act  II 
from  the  "Swan  Lake"  ballet. 


APPENDIX  B 


INTERACTIVE  GRAPHICS  SYSTEM: 

Th.e  graphics  system  used,  for  this  work,  consists  of  a 
cathode  ray  tube  display  and  several  interactive,  devices 
which  are  interfaced  to  a  Univac  1108  computer  with,  a  small 
PDP-8  computer*  The  graphics  system  functions  under  a  swap 
configuration  in  the  1108  „  time-sharing  not.  being-  implemented 
for  the  University  of  Utah."s  1108  at  this  writing.. 

The  display  device  is  an.  IDI  Graphics  Terminal  with 
vector  line,  dashed  line,,  dot  and.  character-generation,  capa¬ 
bilities.  The  interactive  devices  that  have  been  used  in 
this  work  are  a  Syl vania  Tablet  and  stylus  and  a.  Teletype. 
Other  equipment  in  the  graphics  laboratory  includes  a 
Univac  1004  card  processor  and  an  1108  console  monitor-.. 

Software  used  is  described  in.  [201  and.  [2.1t . .  These 
sources  describe  a  general  purpose  set  of  Fortran-callable 
subroutines  for  producing  graphics.  Also  outlined  are  pro¬ 
cedures  for  interacting  by  means  of  the  Teletype  and 
Sylvania  Tablet.  Fortran  V  is  available  under  the  1108 
executive,  the  above  software  packages  being  subsets  of 
Fortran  V. 
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appendix  c 


PROGRAM  LISTING 


is  ASG  DbSUIICCS 

c  ASG  G*S$50?S 
»  XQT  CUR 
I  M  D 

01  FOR  FLEVF> 

LOGICAL  D  F IUI G 
COMMON  OF  BUG. 

0 mu 6  =  .  r  A  L  ?  L  . 

CALL  SYMFM 
CALL  5nn> 

CALL  P) 

r.  m  n 

ON  FOP  S  Y  S  T  F.  *•“ 

SUB ROUT  I  ME  SYSTEM 
CALL  RFI.OAO 
CALL  INPUTM 
CALL  SfTSWP Mil) 

CALL  SFTSlYP  (  M  ,'<00* 

CALL  SrTOUF  (  *4  ,«00» 

CALL  TABABL 
CALL  TMUULMI 
RETURN 
FMO 

ON  FOR  ri 

subroutine  F*  1 

INTEGER  RAGE  I  (  l BOO ) 

INTEGER  X,Y,Z 
NAMEI.IST/N1  /  A  b  (i  R  T  1 
CALI.  STTLST 
R  E  A  D  (  5  t  N  I  ) 

CALL  J'iMTSI  •  N  l  •  fS00) 

I  F  (  F I.  A G  •  G T  •  0 •  I  G  J  TO  MO 
tlagbi  . 

CALL  SFTOr (PAGE  I  ,  I  NOE*  > 

CALI.  PTIIPNI  |  ♦n.o  +  o,  1*0) 

CALL  FI.FINS 
CALL  LCMAR 

CALL  WP| TAT( 3*1  ,flR0, »A  DYNAMIC  BODY  MODEL  •) 
f  At  L  rtRTTATC37S,7S?t*TA*)LE  OF  CONTENTS  *) 

C  A I  L  Sf.HAR 


CALL  “  •  1  I  T  A  T  (  3  6  )  ,  A  '*  P  ,  •  R  I  G  H  T  m  I  P  A  N  G  JL  L  S  •  ) 

CALL  V"-’  !  T  AT  (  6C.C,  ,  A  *•  n  ,  •  ?  •) 

CALL  *  •'  p  I  T  A 1  (  A  5  5  ,  6  5  0  ,  •  3  *  ) 

CALL  ‘-Rl  TAT  <35  1  .HAO, 'MODEL  •) 

CALL  C  f'  I  T  A  T  (  A  6  5 , 4  A  0  ,  •  H  •  ) 

J  5  A  V  E  =  INDEX 
60  TO  60 

HO  INDEX  e  tSAVE 

CALL  RSETDF(PAr;E?  ,  INDEX) 

50  CALL  SfMOr 

CALL  CHR I  NT (  1  ,s70 ) 

CALL  TJMHTU  ,*90) 

60  CALL  IDI  F 
CALL  S  v  A  P 
GO  TO  60 
70  CALL  TTY 
75  CALL  IVTRET 
PO  CALL  SFTOF (PAGE! , INDEX ) 

STOP 

90  CALL  GF.TT  AR  <  X  ,  Y  ,Z  ) 

100  CALL  GF TTAR{ I  ,  J.K) 

I  Ft  K , N  f. -  |  )  GO  TO  100 
irtr.6T.A00)  CALL  P 2 
IF(Y,6T.5l5)  CALL  03 
IF  (Y.GT  •  H  20 )  CALL  PH 
GO  TO  75 
F.NO 


«  TO  AjBORT 


Q  WIPE  HUT  PIC 
ft  it)' HEN  ABORT  TYPED 


0N  FOR  P? 

SUBROUTINE  P2 

INTEGER  X ,Y ,Z ,PAGF? t I «00) 

NAMF.LI6T/N?/A0MRT?,THRH?,!2,ISZ 

CALL  5F.TLST 

Rr.ADt?,  ,02) 

CALL  JMMpst *N2* ,*00, *55) 

1 F  t  F L  A 6 . 0  t . 0  )  6  0  TO  HO 
FLAG*:.  i 

CALL  SFTDF  <  PAGF.2 , 12) 

CALL  COM  A  NO 
CALL  POXES 
CALL  ZL  t  Nf  t  2 ♦ I  ) 

CALL  STICK 

CALL  PTURN  I 2*o  ,  | *0  ,  J  *0 )  , 

CALL  SCHAR 
CALL  FLFINS 

CALL  K R I  TAT { A7 ft  ,  AHO , * RT  HIP  ») 
CALL  005(760, Rift) 

CALL  VI.Cf7I3.729) 

I  5  2  *  1  2 
GO  TO  60 

C  ENTER  HERE  AFTER  FIRST  CALL 


HO 

I2-IS2 

CALL 

RSCTDF ( P AGF2 . 12) 

50 

CALL 

S  F  N  D  (•' 

CALL 

CMfMNT  (  I  ,«70) 

55 

CALL 

TAP  I  NT  t  1  ,  *,90  ) 

60 

CALL 

IOLF 

CALI 

•1,  r  A  P 

5  DRAW  LEG 


ft  TO  ^PORT 
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GP  TO  AO 
70  CALL  TTY 
75  CALL  I*iTRET 
BO  CALL  SFTr>r(PA',E2,I2) 

call  5rM(ir 

STOP 

90  CALL  G  r  T  T  A  0 <  X  ,  Y  ,  Z  ) 

100  CALL  G  f  T  T  A  0 (  !  , J,K  ) 

I  r  f  K  .  N |  »  GP  TO  !00 
iriY.LT.935)  Gn  TO  110 
IFIX.LT.5l2)  CALL  PI 
CALL  P 3 

1)0  ! F ( Y . L  T . a  7 5 )  GO  TO  | 3 0 

i 2* i s; 

CAl.L  G)TCRV(  PAGE  2,  12,2  +  0) 
1 S  HO A  H  12 
C  A  l.  L  C  H  R  1  N  T  (  I  ,  «;  7  0  ) 

CALL  TaB1nT(|,*9(J) 

CALL  IMTPF.T 

130  IF  (  Y  ,LT  .780  )  GO  TO  i  HO 
C  ERASE 

G  0  TO  HO 

I  H  0  !F(Y.LT.72P)  GO  7(1  ISO 
C  SHOT.  CUR VF S 
I  2  a  I  S  H  P  *•’ 

CALL  R$CTDF JPAGF2  .  m 
CALL  SF  OOF 
CALL  I'lTPET 

ISO  IF(Y.LT.A80)  r,0  TO  |  40 
C  RE-USE  CURVES 

CALL  GFTPTS 
CALL  InTRET 

160  IF ( Y.LT.560)  GO  TO  7S 
c  TABLE  OF  CONTENTS 
CALL  Pi 
RETURN 
FMO 


a  >ir>f  out  pic 
rj  WHEN  ABORT  TYPFO 


Q  NOT  page  turn 
©  turn  page 

Q  TURN  PACE 
M  NOT  DRAW  CURVES 

9  DRA'.V  CURVES 


U  NOT  ERASE 


13  NOT  SHOW  CURVES 


a  NOT  TABLE  CONTENTS 


ON  FOR  P 3 

SUBROUTINE  P  3 

INTEGER  X , Y , Z , P AGE  3 (  1  POO ) 

COWMON/US/U (12,9,9) 

N  A  MEL  I  S  T  /  N3  /  A !|  PR  T  3  ,  PLOT  3  ,  !  3  / 1  S3  ,  IJ 

CALL  SFTLST 

READI5.N3) 

CALL  JUUPSI 'N3* ,»P0, $200) 

IF  (  FLAG.GT  t  0.  )  Gfl  TO  HO 
FLAG  ■  I  . 

CALL  SE  TOE ( P  AGE  3 , I  3 ) 

CALL  C  ON  A  M 0 
CALL  PmYFS 
CALL  ZL 1NE ( 3*  I  ) 

CALL  STICK 

CALL  PTUPNI  3*o,  i  *n  ,  |  .r>) 

CALL  SCHAR 
CALL  Firms 

CALL  A  c*  J  T  A  T  (  5  7  A  ,  6  H  0  ,  •  R  T  KNEE  ’) 
CAl.L  P0S(76n,BllR) 
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CALL 

VFC17P3.773) 

CALL 

V  f  C ( 7  5 5  ,  y? 9  ) 

I  S  3  = 

n 

GO  TO 

GO 

ENTER  MERE  AFTER  FIr,ST  C  A 1 

PO 

13  * 

153 

CALL 

p  5  F.  T  D  F  (  P  A  G  E  3  ,  1  3  ) 

50 

CALL 

Sr  NOF 

55 

call 

O'  R  I  N  T  (  |  ,  $  7  0  ) 

CALL 

T  A  5  I  N  T  (  1  ,  «,  9  0  ) 

60 

CALL 

IDLE 

CALL 

5  •  A  P 

GO  TO 

/.0 

70 

call 

tt  y 

75 

CALL 

If  TPfT 

PO 

C  A 1,  L 

s  r  T  0  F  (  P  A  G  F  3  ,  f  3  ) 

CALL 

5  T  0  P 

sr  nof 

90 

call 

G »  T  T  A  A  1  X  ,  Y  ,  Z  ) 

100 

CALL 

GFTTAPII.J.K) 

I  E  t  K  .  N  F  .  -  !  )  G  H  TO  IOO 

if  f  r  .it  .9  35)  on  to  no 
I  F ( X . i  T  .  5  I  2 )  call  p 2 
CALL  P« 

1(0  IF1Y.I  T.P/5)  Gfl  Tn  13  0 
13V I  S3 

CALL  G  F  T  C  R  V  (  P  A  G  E  3  *  1  3  »  3  ♦  0. ) 

15  MOV/  =•  13 

CALL  C  h  P 1  f  J  T (  |  i 5  7  0 ) 

CALL  T  A  0  1  fJ  T  (  1  ,S90) 

CALL  IMTPET 

130  IF (Y.LT.7R0)  GO  TO  190 

C  erase 

GO'  TO  <iO 

IPO  |F (Y.LT.72R)  GO  TO  |S0 
C  SHOW  CORVES 
I  3*1  SHOW 

CALL  R?FTDF(PAF,E3, 13) 

CALL  SFMOf 
CALL  INTRET 

150  lF(Y.LT.ARO)  GO  TO  140 
C  RE-USE  CURVES 

CALL  GFTPTS 
CALL  IHTPET 

160  I  F‘(  Y  » L T  « 560  )  GO  TO  75 
C  TABLE'  OF  CONTENTS 
CALL  Pi 

200  CALL  PLOTor 
GO  TO  55 
END' 

ON  FOR'  P9 

SUBROUTINE  R9 

INTEGER  X , Y , 7 , P AGF9 ( 1  BOO ) 

REAL  LEFT 

COMHON/CPP/PAGEP  f  ISP  ,  l*t 
COMMON/PS/P (  1 7 ,9  ,  I P5 )  , JJ 
COPMOU/M  ISC/NU  ,FY«  STG.'F  YE  TUP 


c  n  m  m  n  r  /  <\  m  g  l  f.  /  t  n  c,  r  1  ?  ,  3  , '  s  > 

NA"tL  1  ST/NS/  A'MIRT  *»  ,  THRUM  ,  !  M  ,  !  5M  ,  fJ  N  ,  F»  ,  FyF.STG  ,  EYE  TUP  .LEFT  ,  IJi 
CALL  SFTLST 
W  F  A  D  f  G ,  M  M  t 

C  ALL  •  JUMPS  «•  NM  •  ,  «5PU  ,  *55  ) 

ir  (flag.gT.o.j  go  to  mo 

FLAG  b  I . 

CALL  S  C  T  D  F  f  P  A  '  T  M  ,  1  M  I 
C  A 1.  L  P  T '  l  R  N  (  M  ♦  r  ,  1  ♦  n  ,  Cl  n  I 
CALL  Lr'MR 
CAl.l.  F  i.  FINS 

CALL  ftp  IT  AT ( 55 , B9A  ,  • SMOrt  MOVE  •) 

CALL  V  »!  T AT ( 5  5 , PM  A  ,  • ROT ATF  X  *) 

CALL  R  l  T  C  T  (  5  S  ,  p  0  n  ,  •  R  0  T  A  T  F  Y  *  ) 

CALL  nJTAT  (55,752  •  ’STOP  •) 

CALL  V.  w  l  T  A  T  (  S  S  ,  S  9  7  ,  •  T  A ,x  L  F  OF  CONTENTS  *) 

CALL  w k  j  T  a  T  f  5  5 , 7  0 M , • R  E - U  S F  •) 

I  S  '*  =  !  M 
GO  TO  SO 

*10  1  4  *  I  s  N 

CALL  PsFTDrjPAGFH.IS) 

SO  CALL  S  E  NO  F 

CALL  C  MR  I  NT  (  1  ,  7  0  ) 

SS  CALL  T  A  0 ! N  T (  |  ,$90) 

60  CALL  I  ^  L  E 
C  A  l  L  S  A  ? 

GO  TO  AO 
70  CALL  TTY 
75  CALL  P-TpfT 
60  CALL  5  r T  0  F ( P  A  G  E  M  ,  !  M  ) 

CALL  SFNOF 
STOP 

90  CALL  GETTAP (  X  ,  Y  ,  7  ) 

100  CALL  GTTTAR ( 1  ,  J  ,K  1 

ITfK.NF_.~I  >  GO  TO  100 
1FCY.LT.9351  GO  TO  110 
CALL  P  3 

110  IFCY.LT.fl75)  GO  Tn  120 
CALL  ACTION 
CALL  INTRET 

120  IF(Y,LT.026>  GO  TO  130 
C  ROTATF  X 

C  UP»V  ALOE  TYPFD  IN  (JU  RADIANS) 

DO  1 25  J» 1 » JJ 
125  ANC,(  1  ,  i  ,  J)«UP 
CALL  INTRET 

130  lr*(  Y.l,  T.  780)  Gn  TO  ISO 
C  ROTATF  Y 

C  LEFTsVALUE  TYPED  (IN  RADIANS) 

on  ns  j=i  ,  jj 

135  A  N  G (  ;  ,  ? , J  )  =  L  E  F  T 
CALL  P'TKET 

ISO  IF(Y.LT«726)  GO  TO  150 
C  STOP 

GO  TO  GO 

150  IFlY.LT.AP,  0)  GO  TO  160 
CALL  GfTPTS 
CAl.  LI  “TRET 
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CALL  P  1 

RETURN 

END 

I3N  FOR  RDATA 

BLOCK  I' AT  A 
REAL  Mu, MX 
INTEGER  P 

C  Cl  M  M  0  N  /  K  1  0  S  /  K  1  [M  1  2 , 3  1 
COMMUM/CONS/flN  (I  ?  ,  3  )  ,  MX  (  1  2 , 3  I 
C  0  M  M  0  /  P  S  /  P  (  1  ?  ,  M  ,  !  **  5  >  ,  J  J 
COMMON /t'LS/UL  (12,3,3) 

COMMON  /M|sC/rjrj,trFST6,rYFnm 
COMMON/IIS/U  (  1  ?  ,  S  ,  H  » 

C  INITIAL  POINT  C  U  0  R  0  I N  A  T  E  VECTORS  (LOCAL  COORDINATES) 

rATA((ncj,<,i),Kxi,M),J=l,‘4)/-Sl?,-S12,-^l2il,  766,700,0,1, 

X  0,-100,0,1,  0,-100  ,0,1/  , 

X  UPCJ,X,ll,Ke|,‘M,J*S,  121/0,  -ino.0,1,  0,-100,0,1,  0,180,0, 

X  -1  10.-2G.0, 1  ,  110, 2S, 0,1,  0, SO, 0,1,  -R  2 , 6  0 , 0  i  1  ,  *(2,60,0,1 

c  ancle  constraints 

X  (  (  MM  (  J  ,K  )  ,  MX  (  J  ,K  1  ,K«  1  ,  3  1  ,  J*3  ,  8  )  /  -2.09880,  .26180,  -1.  309  , 

X  ,  26 1  Bn  .  6 1  Bn , ,52360  ,  0,,  2.0988,  -.261  80,  ,26180  , 

X  -.03,  .03/ 

C  NATURAL  linkage  rotation  MATRICES 

D  A  T  A  (  (  (  U I  (  J  ,  K  ,  L  )  ,  I-  =  1  ,  3  )  ,  K  =  1  ,  3  )  ,  J  «  1  ,  3  1  /  1  .  ,  0  .  ,  0  .  ,  0  .  ,  1  .  ,  0  .  ,  0  ,  0  . 

x  i.,o.,o.|0.,i.,o.,o.,o.|i.,i.(n.|0.,o.,i.,o.,o»,o«,i./ 

OAT  A  (  (  ML  (  J  ,K  ,K  )  ,Kb  1  ,  3  1  ,  J  =  8  ,  1  2  I  /  1  .  ,  1  .  ,  1  ,  ,  1  .  ,  1  .  ,  1  ,  ,  1  .  ,  1  «  ,  1  „  , 

X  ,  1  *  ,  I  .  ,  1  •  ,  1  .  i  i.il.,t<<  1  .  ,  1  •  ,  1  .  ,  1  .  •  1  •  ,  I  *  ,  1  .  ,  1  .  ,  1  .  / 

c  tree  structure  or  model 

D  A  T  A  (  (  k  1  D  ( »> ,  K  )  ,Ksl  ,3j  ,  .1  e  1  ,  121/0,0,2,1  , 0, 3, 2, S, 8, 3, 0,0, 2, 7,6,5 
X  7,0,8,7,9,0,7,10,0,7,0,11,10,12,0,10,0,0/ 

c  m i sc ,  constants:  nn=«  p t c  thinned  to,  nh  ,le.  25, 
c  EYESTCi  AND  EYETuBs  ARS.  DISTANCE  TO  STAGE  AND  DISPLAY  TUBE 
DATA  N  u  ,  E  Y  E  S  T  G  ,  f.  T  E  T  U  B  /  2  r»  ,  2  0  ,  ,  1  •  M  /  1J  (  1  ,  8  ,  8  )  /  1  ,  / 

END 

ON  TOR  G E  T  C  P  V 

SUBROUTINE  GETCRVCOE,  IMOF.  X,  IP) 

C  THIS  READS  3  CURVES  IN  ANY  ORDER  AND  STORES  X'S  IN  5  A VE X ( JO  I  NT , BP 
C  AFTER  HITTJNG  ORA -CURVES  LIGHT  BUTTON, 

C  ONE  OP  MORE  CURVES  MAY  BE  P  E  0  R  A N  » 

C  AFTER  RE-OWAfl  1  f|Ci  ,H  j  T  KEEP-CURVES  LIGHT  RUTTOO, 

logical  debug 

INTEGE  R  OF ( 1 800) , INDEX , X , Y ,7  ,  TFMP 
COMMON  DFO'IG 

COMMON /CSX X  5  A  V  E  X ( 3,1001,SAVFY(3,lfl01,IS(3l 
N A *•’ E L  I  ST/NC  /  ABORT 
CALI  SETL5T 
R  E  A  0  (  5  ,  N  C  ) 

CALL  JUMPS ( • NC» ,$1001 
CALL  C H R 1 N T ( |  , $ 9 0 ) 

IF ( IP. CO,  IPSAVE  )  GO  rn  ? 

IPSAVEr IP 
OP  1  K* 1 , 3 

1  I S ( K )*0 

2  CM  M  TRY  BOX  It  I 

5  CALL  TABlNT(l,s201 
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10  C  A  l.  L.  1  ri  L  E 
CALL  S'.- AO 
GO  TO  10 

20  CALL  fir  TTaB ( X ,y , / )  y  GET  A  POINT 

I r ( y . i  t . p 2 n  )  go  to  ?  \ 

I  F  1  {  Y  .  G  T  .  ft  7  S  1  .  n  R  .  {  x  •  Cl  T  .  5  I  ?  )  )  GO  TO  2  1 
C  LIGHT  PUTT’  N  FOR  >  >  Cf'  CURVES 
CALL  GT  T  ANG  <  I  <>♦  1  ) 

C  IF  EXPAND  I  MG  PROGRAM,  CHANGE  FOLLOWING  INTEGER. 

I  F  1  I  P.f.O  .  3  )  CALL  Gf  TPTS 
2  2  CALL  GFTTAB(X,Y,Z) 

I  F  1  7.  .  N L  .  -  1  )  GO  TO  2  2 
GO  TO  SO 

21  TF.MP=<k~1M3S6 
jt=n 

I  F  (  X  .LT  •  (  G'l  +  TFMP  )  .OR  ,  X  .r,T  .  (  ?5A*TLMP  J  )  GO  TO  *10 
GO  TO  2 o  9  RIGHT  BOX 

27  CALL  GFTTAB ( X , V , 7 )  «  GET  A  POINT 

2 R  t  F  (  7.  .  F  <'  .  -  1  )  GO  TO  3  0 

I F  (  Y  .  1.  T  .  o  0  .OR,  Y  .  G  T  .  1 A  2  )  GO  TO  27  (3  OUT  OF  RAnGF 

IFTX.LT.  CA*1  +  T1MP)  .DR.X.GT.  (75/.  +  TEMPJ)  GO  TO  27 
J  -  J  «■  I 

SAVEX  1  *'  ,  J»  =X 
SAVE  Y  (  <”  ,  J  1  »  Y 
GO  TO  2 7 

C  DISPLAY  DOTTED  LP‘F 
30  J  S  (  K  )  B.l 
I  s  I  NDf  < 

CALL  RSFTDF  (OF  ,  I  ) 

DO  3  5  Ys  )  ,  3 
,1b  I  S  (  Y  ) 

IF (J.Ffi.O)  GO  TO  35 
DO  3 *1  X »  1  t  J 

3*1  CALL  DriT(lFlX(SAVF.XCY,X)),iriXlSAVEY(YfX)}) 

35  C  TINT  IN' IF 

call  s *: m f ) f 

CALL  IMTRET 

*io  KnK*i  o  try  another  n ci x 

IF  (  K  •  F r'  •  *1  )  K  =  ! 

CALL  TnTRET  0  RE-ENTER  AT  $20 

C  all  THROUGH 

50  JNDEXbJ  «- 

IF(DEHljG)  CALL  T  YPOUT  (  •  GOTCRV  •) 

RETURN 
90  CALL  TTY 

CALL  If. TRET 
C  EMERGENCY  EXIT 
100  STOP 
END 

ON  FOR  C,r  T  ANG 

SUBROUTINE  GETANG 1  I  PAGE ) 

C  FROM  ARRAY  XN,YN  THIS  COMPUTFS  THREE  ANGLE  ARRAYS 
C  THIS  SUBROUTINE  USES  ANGLE  CONSTRAINTS 
LOGICAL  DEBUG 
REAL  M'!»MX|M,LX 
COMMON  debug 

COMMON/PS/P(  I  2  ,*1  ,  I  *15  I  ,  JJ 


48 


c 


COMMON /  TNS/XMO.INM.YU 
COMMON/ ANGLE/ A N6 ( 

C  0  M  M  0  N  /  C  0  N  S  /  M  M  {  12,3)  » M X  < 12.3) 


CALL  l’-'TRP? 
00  20  *•=  1  ,3 
pr  20  J* i  .  J J 


lol  13  0  X  B  =  AX1S  ABOUT  WHICH  ROTATE 
»  I M DIVIDUAL  DATUM 


LX*AM+(K'-|  J  *  3SA 
RXt!2SF>-*.  (  <-  1  )  •  X&A 

M  =  L  X  +  <  A  R  S  C  M  N  (  I  P  A  G  E  » K  1  )*)92)/(ABS(MN(  I  P  A  G  F.  *  K  )  )+MX(  I  P  A  6  E  »  K  )  ) 


1  T  (  X N  l  y  .  J )  •  G E  •  M  )  G U  TO  10 
NEGATIVE  A.iGLE 

A  MG  (  IP  AGE  *K  ,  J  ).  «  C  M-XNIK  %  J.  >  >  *MN  t  |PA5E|K)/<M-LX) 


GO  TO  ?0 

C  POSITIVE  ANGLE 

JO  ANG( IPAGE.K, J)«( <XN(<,J>-M>.MX(  IPAGF.K)  » / «  RX-M  > 
?0  CONTINUE 

IE  I  DEBUG  )  CALL  T YPOUT ( • 60T ANG  •) 

RETURN 

END 


CN  FOR  GfTPTR 

C  UPDATES  MATRICES  AMD  FINDS  ARRAY  OT  DISPLAY  VECTORS. P. 
C  CALLS  MATMUL  ,  VECMAT  AND  ROT.  WHICH  CALLS  M3X3M. 

C  DRIVEN  DY  TMf  TABLE  ENTERED  IN  PLOCK  OATA. 


SUBROUTINE  GETPTS 
INTEGER  P 
LOGICAL  DEPUG 
COMMON  DEBUG 

C0MM0N/PS/PCI?,Ril‘*5)»JJ 
COMMON  /US  /U  (  1  7  ,H  , ) 

COMMON/E IDS/K I D I  1 ? i 3  I 
COMMON /M I  SC/ Mu .EYESTG.EYETUB 
CALL  TvPOUT ( M 1 3HEYE  TO  STAGE* 
X  FS.2.  ?H  )  •  .EYESTCuEYrtUR ) 

EToEYETUB* I  000. 

P( l .3.1 1=“EYESTG*70» 
c  TABLE-DRIVEN  POINT  DETERMINATION 
DO  30  .1*7  .  U  J 


ES.O.SM  FT#  I6HEYE  TO  DISPLAY  » 


ID  NUMBER  OE  ’FRAMES* 


NODE*  1 

C  STORE  ROTATION  PART  OE  MATRIX  IN  UPDATE  MATRIX 
18  CALL  POT  I  NODE  » .J ) 

HOAD**  j  D  (  MODE  ,  1  )  •  ■  . 

I F ( NODE • EQ • I )  GO  TO  22 
C  UPDATE  ROTATION  PORTION  OF  MATRIX 
CALL  MATMUL (NODE, NO  AO) 

C  FIX  remainder  DE  matrix 
DO  20  !  =1 
U  (  NODE  ,L  i*H  *0. 

20  II  (  NODE  |L  >  »P  (  NDAD  ,  L  .  I » 

C  GET  POINT 

22  CALL  VECMAT ( NODE i J1  . 

C  HERE  COULD  BE  DETERMINED  COORDINATES  OF  ALL  STRUCTURES 
C  THAT  ARE  FIXED  IN  THIS  COORDINATE  SYSTEM. 

C  TEST  FOR  SON 

IF (KID (NODE. 3) .EQ.OJ  GO  TO  ?A 
NODE»K 1 D ( NODE , 3 ) 

GO  TO  1 R  ' 

;  TEST  FOR  BROTHER  ■ 
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2ft  IF  <K  10  (  NOpF.  ,2  >  .F^.O  )  GO  TO  2ft 
MfjnpeK  I  n  (  NODE  ,  p  ) 

60  TO  1  n 

C  TEST  FOR  f  ATMER 

20  1  F  (  K  f  0  (  MODE  ,  1  >  .  E'.-i.O  1  60  TO  30 

C  back  UP 

n  n  o  e  *  *  i  o  c  n  n  o  f  ,  i  j 

GO  TO  7ft 

30  ClJNTIU'iF 

C  PFR5RECT1VF  TRANSFORMATION 
DO  32  J=?,JJ 
DO  32  1  P  =  2  .  ! 2 
DO  31  1=1,2 

31  P(|P,L,.l>=FLnAT(PnP,l.  ,J))*LT/FL0AT(-P(lPt3tJ))*5|2, 

32  CONTINUE 

I F ( 0  E  R  n  6  j  CALL  T V P OUT < • GO T P T 5  •) 

RETURN 

END 

ON  FOR  MATMUl 

5  U  R  R  0  U  T 1 N  E  MATMUL (NR.NA ) 

C  UPDATER  ROTATION  PORTION  OF  MATRIX 

C  MULTIPLIES  MATRIX  ll  ( NP )  X  n(NA)  AND  STPRES  IN  U(N0) 

DIMENSION  A ( 3 , 3 1 
COMMON /  US  / U  (  1  7  ,  ft  ,  <1  ) 

DO  I  .1*1  ,3 
00  I  K* I ,3 

1  A ( J  ,  K 1 =0 • 

DO  2  J=l  ,3 
00  2  K* I , 3 
on  2  L  =  1  ,3 

2  A(J,K)=A(J,<)  ♦  U  (  N  B  ,J  ,  L  )  •  U  (  N  A  ,  L  ,  K  ) 

DO  3  .1=1,3 

00  3  Ks I ,3 

3  U(NB,.I,K)  =  A(J,<) 

RETURN 

END 

QN  FOR  VFCMAT 

SUBROUTINE  VFCMAT IP, | } 

C  MULTIPLIES  INITIAL  ROW  VECTOR  X  MATRIX  M  AMD  STORES  IN  VECTOR  P( 
INTEGER  P 

COMMON/US/UU  ?  ,R  ,Rj 
COMHON/PS/Pl  12, ft,  MS)  ,JJ 
DIMENSION  V(ft) 

DO  I  L=1,M 

1  V(L1=0. 

DO  3  K  «  1  ,  ft 

on  2  L « I  ,ft 

2  V(K)=V(K)+ELOAT(P(M,L,I  )  )*U(M,l.  ,K) 

3  P (M , K , I )  =  1 F I  X (V  IK  14.00051 
RFTURN 

ENn 

BN  FOR  POT 

SUBROUTINE  RfTT  (  JOINT,  J1 
C  SFTS  UP  AND  MULTIPLIES  ROTATION  MATPJCE5. 

C  PLACFS  PRODUCT  IN  UPDATE  MATRIX,. 
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LOG  1C  At  DEBUG 
common  DEBUG 

DIMENSION  D ( 3 , 3 1  •  E  (  3 , 3  1  ,  F(3,3) 

COMMON/ULS/UL <  t  2 , 3  *  3  1 
COMMON /UG/U  (  1  2 ,H  ,'U 
COMMON/  ANGLE/  ANG  (  1  ?  ,  3 , 1  ‘15  1 

c  SET  UP  ROTATION  MATRIX  FOR  NATURAL  LINKAGE  SYSTEM. 
DO  1  L  *»  l  i  3 
DO  1  K  *  1  ,  3 

1  F  (  L  .K  1  =UL  (  JO  1  ft  T  ,L  ,K  ) 

C  SET  UP  X  ROTATION 

2  D(l,l)=l. 

0(1  ,  2 1  rO. 

0(  1  ,31=0. 

0(2,1 1=0. 

DC2,2)»C HSIANG (JOINT, I  ,J)  1 
0(2,3)=SlN(ANG(jniNT,l  ,J)  ) 

0(3,11=0. 

0(3,2)  =-0 (2,3) 

0(3,31=0(2,2) 

C  MULTIPLY 

CALL  M  3  X  3  M ( F  ,  0  j 
C  SET  UP  Y  ROTATION 

E ( I  , I  1  =  cns( ANG (JOIN T,2,J)  ) 

E( 1  ,21=0. 

E (1  ,3 1 =-S 1 N ( AMG( JO  1  NT ,2  ,  J) 1 
F ( 2  ,  1  1  =0. 
r  ( ? ,  2 1  =  t .  . 

F (?,3)=o, 

f  (  3  ,  i  )  =  ~  r.  (  i  ,  .1 ) 

F.  (  3  ,  2  1  *  0  .  " 

F(3,3)=F(1,I) 

C  MULTIPLY 

CALL  M  3  v  3  M ( D  ,  F  ) 

C  SET  UP  7.  ROTATION 

0(1,11=  COS ( ANG ( JO  1  NT  ,  3 , J)  1 
0(1,21*  S IN ( AUG ( JO  1  NT ,3  ,  J)  1 
0(1,31=  0. 

0(2,11=  -F (  1  ,  2  ) 

D(?,2)«FII»1) 

0(2,31*0. 

0(3,11*0. 

0(3,21*0. 

0(3,31*1. 

C  MULTIPLY  ROTATION  MATRICES 
CALL  M3X3M(F.,0) 

C  STORE  IN  UFOATE  MATRIX 
[)0  JO  1*1,3 
00  10  K  *  1  , 3 

10  U( JOINT ,K ,L )  =  0 ( K  ,L 1 
I  2  RETURN 
ENO 

ON  FOR  M 3 x 3 M 

SUBROUTINE  M  3  X  3  M  (  I*  ,  A  ) 

C  MULTIPLIES  MATRIX  B  X  A  ANO  STORFS  IN  A 
DIMENSION  A(3,3),0 !3,3),C(3,3) 

00  I  J= 1 , 3 
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DO  1  Kn 1 , 3 

1  c  <  j  ,tc  i  =o. 
no  7.  J=1  ,3 
DO  7  K  = I , 3 
on  7.  L  s  |  ,  .1 

7.  C  (  J  ,K  )  rC  (  J  ,K  )  +f>  (  J  ,  L)  *  A  (  L  ,  K  ) 

00  3  Js|,5 
on  3  r.r.\  ,3 
3  A(J,K)~C(J»K) 

RETURN 

END 

I3N  FOR  ACT!n*. 

SURRnilT  1  ME  ACTION 
C  DISPLAYS  NOTION 

C  DRIVEN  B Y  THE  TABLE  ENTERED  I M  BLOCK  DATA, 
INTER  I  PAGER! |ROP)tp 
C  0  M  M  fl  N  /  c  r  H  /  P  A  r,  f  M  ,  I  s  h  ,  m 
c  n  p  m  n  n  /  n  s  /  p  n  ? , «» ,  i  r  5 ) ,  j  j 

COMMON/*"  I  DS/K  I  n  (  1  7 ,3  ) 

CALL  R  5  F  T  f)  F  (  P  A  6  L  ‘I  ,  I  H  ) 

on  3  J=2,JJ 

IlrlSB 
NF  r  p 

CALL  PPSIPINF , J , J) ,P(MF ,2, J) ) 

10  NSeK 1 P ( NF , 3 ) 

IF (NS.rO. 0)  60  TO  IS 
12  CALL  vr.c  (  P  (NS  ,  |  ,  J  )  ,P<  MS  ,  ?  ,  J  )  ) 

N  F  r  N  5 

go  Tn  in 

15  NS«K 1 0 ( MF , 2 ) 

IF (  N  S  •  F  Q  •  0  )  GO  TO  20 
NFcK 1 0 ( NS  ,  1  ) 

CALL  POS (P(NF, | , JJ ,P(NF,2 , J!| 
cn  to  i? 

20  IF (K  ID(*T  , 1 ) . FM.O)  GO  TO  25 
OF eK 10 (NF  ,  1  ) 

60  TCI  15 
25  CONTINUE 

CALL  v  r;  C  (  p  (  1  1  ,  1  ,  J)  ,P(  l  1  ,2  ,  J)  ) 

3  CALL  Sr.MDF 
RETURN 
END 

ON  FOR  BOXES 

SUBROUTINE  fJOXFS 
CALL  FLE1NS 

tin  i  *=i  ,3 

CALL  55  BOXES 

CALI.  Vff  ,  (|Q| 

call  Vi  C(7Sfi-.l,*i*2h 
CALL 

CALL  VF  C  I  J*  if] 

CJLL  PnS (55VJ.1T5 J  «  ARROWS 

CALL  Vf r I  55- J, | 

CALI.  VTC  l“T  *  J  ,  J  1  fl  1 
CALL  t**S I  55- J , ) JO J 
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CALL  vec  C  A3«-J  ,3)  8  > 

CALL  SC*<AM 

CALL  F8  I  TAT  I A8  +  J, AO, • min  •)  «  LEGTND 

?  CALL  *»tTAT(23?«.J,AO,»MAX  ’) 

CALL  LCMAR 

CALL  W.MTAT(  177, 870  , *X  *  ) 

CALL  *i?!TAT(5?8,H?0,*Y  *) 

CALL  *'  P  l  T  AT  (PM  ,870  ,  •  7  *  ) 

RETURN 
f  NO 

IJN  rOH  5TIC< 

SUPROUT  IMF.  STICK 
C  FIGURE 

CALL  FI.EIN5 
CALL  P  N  S  (  7  A  fl  ,  9  0  1  ) 

CALL  VTC  I  78H  ,<>  30  ) 

CALL  V  f  C (780,930) 

CALL  vrf(7A8,9ni) 

CALL  OPT  I  7 A  I  ,9?3 I 
CAIL  0OT(775,°J>3) 

CALL  P  0  S  (  7  A  5 , 9  I  8  ) 

C ALL  Vr  C  I  77  I  ,9  I  8  ) 

CALL  P'St7A0,901 ) 

CALL  VrC 1760,708 ) 

C  AXES 

CALL  PflSIP?fi,57A)  5»  X 

CALL  vrC  I  7Afi  ,S7A  I 

CALL  vrCC7An,A3A)  *>’  T 

CALL  P M f»  (  7 A 8  , r? 7 6  ) 

CALL  vrC  «  750 ,588  )  ®  Z 

CALL  Sf.HAR 

CALL  vM?  I  TAT  (W. 3  0,57?  ,  •  X  *) 

CALL  «'MTAT|7^,A'U,'*  •  ) 

CALL  */R!TATl7?f  ,**19,*Z  •) 

RETURN 

f.np 


ON  FOR  7LINE 

C  ENTERS  DAS*'F0  LINE  OF  DISPLAY  POXES  IN  DISPLAY  FILE, 
SUPR0UT1NL  ZLIME(V) 

CALL  !  T  AT  c  35  I  ,55  0  ,  •  R  I  GMT  KNFF  ANGLES  •) 

RFAL  r-: ,  rX 

CONHON/CNNS/MMI  |?,3)i’»'(II7,M 
CALL  Fl.EfNS 

X«»llS  (  mm  C  K  ,  II  I  /  (  A  ns  <  MNI  K  ,  I  )  )  ♦  **  X  (  t  ,  1  )  ) 

Y*AOS(f  N(k ,  2  ) )/(APS(MNIK ,2) (♦NX  IK  ,?!  ) 

ZaAOSI  -  ‘A1  (  K  ,  3  I  )  /  (  A  P  S  (  M  N  (  K  ,  3  )  ) +HX ( K  ,31  ) 

J*X* 1 9?* 68 
CALL  P'1S»J,P0) 

CALL  0ASIMJ.MA2I 
J«Y*19?*A8*35A 
CALL  PPSIJ.OO) 

CALL  PASH(.‘,8A?) 

«A8«2o 5  A 

call  rr s i j,poj 

CALL  D  A  S  H  (  J  ,  8  A  2  ) 

RtTUR*i 


53 


END 


fiON  FOR  Chmanp 

5UBR0UT  I  *'F.  CO  ■‘AND 

C  ENTERS  COMMAND  QUA  f)  RA  NT  IN  DISPLAY  FILE 
CALL  LC"AR 
CALL  F|.  F  I  ns 

C»LL  *  M  I  T  AT  I  54.,  *  ITCtAft  £imv£9  tj 

CALL  "^MAl^^hHn.tfFrr  fLIPVft;  i} 

CML  C^ITAT  |  S<i ,  nrj"t  If  If  ASF  >> 

CALL  **•  IMTi*,*  t  2S  ?  *  •  \ 

CALL  *MTATiSS|jn*(*Rr-IKE  -J 

CALL  Hi  t  TAT  ISS.SWJ  *  *UMt,F  Of  tOMTFNlS  <| 

M£TLJR«i 
t  Nfl 

QN  FOR  PTURM 

SUBROUTINE  RTIIRN(NIIMP,LEFTA,RIGHTA) 

C  ENTERS  PAGr  NUMBER  AND  TURN  IN  DISPLAY  FI|E 
I  N  T  E  fi  F  9  R  I  OH  T  A 
IF(LFFTA.EQ.O)  C,{)  TL1  ? 

CALL  F|.  F  I  NS 

CALL  Pns ( M«fl  ,9AA  )  B  LEFT  ARROW 

call  vrc<97?,96M 

CALL  vrC(M77,?70I 
CALL  PrS ( 97? ,V6A ) 

CALL  VFC{977,96?) 

2  IFIRIGmTa',  EQ.nj  GO  TO  3 
CALL  FLFJNS 

CALL  P(1S(S.16,9AA»  0  RIGHT  ARROW 

CALL  VrC(SS2t9<,6) 

CALL  VFC  I  S*t7 ,96?  I 
CALL  Pnsisr»?,9A6) 

CALL  VrClS'17,9701 

3  CALL  SCUAR 

call  fleins 

CALL  YR I TA T { ,  V6  I  ,  •  TURN  ») 
call  lchar 

CALL  «P!TATf<l96l9ft2l>(?Hp(,f?||H  )*,NIIMP  ) 

RETURN 

END 

ON  FOR  EVEN 

SUBROUTINE  EVEN 

C  FINOS  ARRAY  X,  EQUIDISTANT  BY  SCOPE  y»s 
C  Incoming  apray  savex.save*,  each  is  long 

C  OUTGOING  ApPAT  X,Y,  EACH  "N  LONG 
C  SAVEY  IS  SE  ALEO  TO  RANGE  BETA,  l-NN, 

C  THEN  thinned  TO  VALUES  CORRESPONDING  Tn  SCALED 
C  VALUES  CLOSEST  TO  THESE  INTEGERS. 
logical  DEBUG 
common  nrnUG 

COMMON /CSX /SAVE* ( 3 , 1  BO ) »5AVEY( 3,180)  ,  15(3) 
COMMON/ X S/ X ( 3,75)  ,  Y(3,?S) 

COMMON /M ISC /NM , EYES TG,  EYE  TUB 

on  ?n  npoxb t ,  3 

VfMDOX ,J  IcSAVEYI  NB  0  X  ,  I ) 

X (NBOX , I ) »SAVE X (  NBO  X , |  ) 
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KKn  I 

I N« I S ( rnox l 
no  10  ¥  =? ,  f  n 

Y  M  0  R  M  a  (  R  A  V  E  Y  (  NHOX,<  I  -  «>  O  .  ).NN/3fn 

1T(  ynopm.lt.fi  oat  hoc  »  >  r,o  Tn  in' 

KKnKK+1 

Y  (MDOX  ,KK  )=  5  A  V  E  Y  (  NBnx.K) 

X  f  UFI0X  ,FK  )  *SA  VE  X  f  NBHV.K) 

10  CONTINUE 
20  CONTINUE 

IF  I  DEBUG  )  CALL  T  Y  POll  T  (  •  r  VE  NF  n  .. 
PETURN  ' 

FMD 


S3  R'JNS  FROM 


I»NfJ,  Mp  y 


I  N  T  E  c 


them  filled  |h  A  G  A  IN, 


ON  FOR  INTRP2 

sun r n u time  i  n  t  r  p  2 

C  POINTS  E  YE  ML  Y  THImmeD, 

LOGICAL  DEPUG 
COMMON  D f D 0 g 

COMMON/ X  MS/ XN  I  3  ,  I  RAI  ,  YUJ  3  ,  m6, 

COMMON/ X 5/ X( 3, 25) ,Yf3,?5| 

COMMON/GAUS/H  <  3  ,  M  I  ,  NO  ,  A  (  «|  ) 

C0MM0N/LS7/N.MP0X 
COMMON/PS/PI  12, M, IMS)  , jj 
.COMMON /M (SC/ NN ,E YES TG ,  E YET Uo 
IFfOEPllG)  CALL  TYpOUT(»OFGlN  INT2 
CALL  E  V  F  M 

FIT  EACH  M  PTS  TO  LST  SODA *F  PARANOIA  AND 
INTERPOLATF  5  PTS  BETWEEN  2  AND  3.  ° 

FIRST  AND  LAST  INTFRVALS  APp  TRFATrn  . 

n«*i  ‘  0  RV  linear  interpolation, 

.  IaI  Wn.  DTP  IIPr*A»  _ 


) 


c 

c 

c 


0  Nn •  PTs  USED  TO  FIT  parabola 


C  JJ 


9 

V 


<»• 

n  o  x 


jj= 

« 


NO,  FRAMES  **» 


w  NATCH  npoes 


20 


20 


30 


N«M 

NNM  a  mN-i 

NNH2«N\M- ! 
computed 

J.ln  {  Nfl-  I  )  *A*  t  x 

00  30  jr.  |  ,  3 
N0n3 

On  20  I  Jr  |  ,  NNM  ,  NIJM2 
I JPn I J*| 

00  20  KnIJtIJP 
I  Ka  (  K«*  t  )  «  A  ♦  I 
YN(  J  ,  I  *  )  *Y  (  J  ,(f  ) 

XN(J,lF)»X(J,K) 

00  30  x*l | NNM , MNM2 
DO  30  L=?|A 
Y?«Y( J,F*| I 
YI  «  Y ( J ,  K ) 

I  F  »  AOS  (  Y2«  Y  |  |  ,  r,T  •  I  2  ,E- 1  3  )  I  GO  Tn  ?B 
IFIDEDOG)  CALL  T Y rOUT ( • Y2-Y | =o  i, 

Y2«Y2*2»L-I3 
I  I  * ( K« | I • A  *L 

YN(J,  I  |  ) «» Y  I  ♦  (  L  ••  I  •  )»(V2“Y)  )/A, 

XM|JtIlIa(YN(JtII)-r|,.IX(J  K  ril/f¥,  VI 

00  MO  N  P  n  x  a  l  ,3  J,F  )  )/(Y2-Yi  )  ♦  x  <  J  ,  If  | 

DO  MO  K  * 3 , NNM 

call  lstsqim 

I  I  a ( K- 1  ) • A  ♦ I 
YMfHOOX , I  I )>Y(NBOX ,K I 


®  EII.L  IN  5  PTs 
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JJIJilJ..  JIW.  ILUI-U 


XN(NROX I  I  I  )  =  K I NBOX  ,K | 

on  mo  .ib  >  ,5 

L«=! 

YN(Nnnv,L)«YCiPOX,K-|  >  ♦  J*  (  Y(NPOX,K»-Y(NnOX,K-| 

TEMRbY/'M  NRnx  ,L  ) 

oo  xni  nnnv  tn  nA  m  i+a(2)«tfmp+a(3>*temp»tfmp 

L  a  fj  N  M  *  /,  ♦  1 

90  RETURN 
END 

QN  POP  LSTSG; 

SUOPOUT I NE  LSTSC(K) 

C  SCHAlJM  ST  A T  l  ST  I CS , SP I EGEL »p • 22 1 
C  GET  CONSTANTS  FOR  X  =  A  1  ♦  A  2  •  Y  +  A  3  *  Y  •  *  J> 

C  N  r,  NUMBER  OF  PniMTS  USED  IN  FITTING  APOVE  PARABOLA 
C  X  IS  OFPENoEMT  VARIABLE.  K  |S  3RD  PT.  OF  PARABOLA. 

logical  debug 
common  debug 

COMMON  LST/N.moox 
COMMON /G A US/O  r  3 , M )  . NO , A ( R ) 

COMMON  'XS/X  (  3  ,  ?*>)  •  Y  (  3 ,25  > 

C  Find  COEFFICIENTS  OF  »>OPriAI.  EQUATIONS 
DO  10  jt  |  ,3 
00  10  l,.r  I  ,*< 

10  R(J,L)  r  Q. 

0(1  •  I  I  r N 
DO  20  Jr  I  ,N 
XPr.XI  WiDX  ,K.  J-3  ) 

YP  n  YfMPOX ,K*J-3) 
n  (  I  ,  2  )  n  n  (  )  ,  7  )  4.  Y  P 
0(  I  .MlrOl  |  .MI+XP 
B(2i<|)bR(2,<|)4.VPbXP 
YP«YP**vP 

n<  l  ,3)«:n(  i , 3 )  *  yp 

0  (  3  ,  *»  )  rB  (  3  ,  B  )  +YP«  XP 
YPbYPwy {  N  0  0  X , y 4 J-3 I 

n  ( 2 , 3 1  o  c  r ,  3 1  ♦  y  p 

YPnYP*  Y  (  MO  OX  ,*'♦  J-3  ) 

20  0(3,3)rO(3»3) *vp 

B(2,2)rO( i ,3) 

B  (  3  ,  1  )  r,  0  (  I  ,  3  I 

0(3,21*0 <2, 31  , 

0(2, I  )«0( I ,2) 

C  find  A*s  nr  gauss  fl IMINAT ion,  MACK-SUBSTITUTION. 

CALL  GAUSS 

return 

END 

QN  FOR  GAUSS 

SUBROUTINE  GAUSS 
C  WFFG  AND  ‘■’F  FD  ,  P.  |  OP 
C  SOLVES  SYSTFM  PF  F.nUATIOfJS  BY 
C  GAUSS  ELIMINATION  AND  H  A  CK -SU*»ST  I  TUT  |  ON  • 

LOGICAL  0 L 0 U G 

common  ornuc, 

COMMON /G  AOS /B  (  3  ,  *n  ,  NT  ,  A  <  *4  ) 

NPbNP*! 

N  M  a  M  p  -  I 


)  )  /b  . 


\ 
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c  REDUCE  Tn  T  R  I  ANGULAR  F  f?  R  M 

DO  10  Jr.  1  ,  MM 

JPa J* 1 

0  0  |0  r=JM,MH 

p*n i k  ,  1 1 /ft i J , j ) 

on  in  i r,j , ‘j p 

in  nor  ,i  isnif.i  i-pmiij,  1 1 
c  solve  $  y  r.  t  r  ,f  p,  y  n  Af  k  -suhst  t  tmt  i  nn 
00  20  v  =  NP  ,1,-1 
?0  AIKInfl. 

00  3  0  I'nfin,!  ,-| 

K  r*  =  K  *  1 

00  25  .lr"P,KP,-| 

25  A  (  K  1  r  A  (  K  1  -0  (  <  ,  J  )  *  A  (  J  J 
30  MKlsl»irtNP)  +  A|K)l/!'K,YI 
nrTupr! 
ton 

13  I  FOR  SOCK1 

SUP  ROUT  IMF  sun." 

C  THIS  DISPLAYS  T1TLFS  USED  IN  HALING  TMF  M  p  V  I  F. . 

INTEGER  OF  I  I  ROp  ) 

N  A  M  F.  L  I  S  T  /  5  H  /  1  ,  G  2  ,  S  3 , 6  «  ,  G  5  ,  G  A  ,  G  7  ,  G  P  ,  G  9  ,  G  1  0 , G  1  1 
CALL  SrTLST 
RF  A  0 ( S  ,  SH I 

CALI.  JUMP'!  (  ,  s.h  ,  ,,,]  l?>/»,«:3,1,‘4,-«,S,*A,S7li:R,59,S|0,$lI  ) 

CALL  f r Tr:r  (  OF  ,  I  } 

CALL  CMP  I  UT(  J  ,  <;3n» 

CALL  t.  r  1 1  A  i? 

70  CALL  IPl.r 
CALL  S • A  P 
60  TO  ?0 
30  CALL  TTY 

CALL  I’  i  T  p  F  t 

1  I"0 

CALL  Y<r  I  TAT  (  ,39?  ,  A«M  ,  *OYMAM  I  c  STUDIFS  M 
CALL  m.|i  j  T  A  T  (  9  <IP  ,  AIM  ,  •  FOR  •) 

CALL  op  f  T  AT  (  ?o*  f  f ,  CPMPllTrR- MDEO  CHOREOGRAPHY  •) 

CALL  i-.  P  !  T  AT  (  M(i,i  ,  hSO  ,  *  C  A  PPL  '."I  THROW  ») 

CALL  apIT  A  T  (  27?  ,  MOO,  •  Of  r  ARTHEL’T  IH7  COMPLITFR  SCIENCE  •) 

CALL  v*  R  I  TAT  |  3AH  ,3G0  , »  0**1  VFPS  I  T-Y  f.p  UTAH  'I 
CALL  Y'P  I  T  AT  I  ROR  ,  3PO  ,  t  pf«RUARY  L*>70  •) 

C  A  |  L  S  r  H  0  F 
C  A 1  L  I  OTf’I'T 

2  I«0 

CALL  V  RJ  T  AT  (  ,  c,nM  ,  •  THF  END  M 

CALL  S'7  "OF 
C  A I  L  If.'TPET 

3  laP 

CA1L  I  T  at  (?<!/,,  952  ,•  THE  Mn  V  r  a.  r  M  T  is  SrOTM  IN  2  5  FRAMES.  *1 
CALL  Sr’iOF 
C  A I  L  I  "TPFT 
1  IbO 

CALI-  i*  R  ?  T  A  T  (  3  5  ?  ,  5  0  0  ,  •  T  H  F  SAr'F  MOV  Ff'F.  MT  IS  M 
CALL  V  l  • 1  T  A  T  (  3  c.  ?  ,  **'  S  ?  ,  ’  S  l 1 0  '  f-l  IN  20  FRAMES.  M 
CALL  SFORF 
CAM.  H'TM'T 
5  \r.n 


/ 
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CALL  '.VP  \  T  AT  (  V.  r  ,  r,nn  ,  •  THf-  s  A  M  F  MOV  EM  NT  IS  •) 

CALL  ftPITA":  (?s?1MS?,'Scn..iri  in  i  a  FRAMES.  •) 

CA!  L  STUFF 
CALL  PiT^ET 
A  Inn 

CALL  p  J  C  I  I  •  *1  ,  V) ,  0  F  ,  I) 

CALL  IMTI’CT 
7  I  rfl 

call  p  i  c  ( i  •  *4 ,  ?  n » n  f  1 1  j 

CALL  I'.  7  PE  T 
ft  Inf) 

C  A I.  L  P  |  f  (  3  .  ,  3  0  ,  D  F  ,  1) 

CALL  p'TKFT 
9  I  n  0 

CALL  ft  K  I  TAT  (  3  9  o  ,  S  ft  H  «  •CAMERA  ROTATED  •) 

CALL  5  r  M  D  F 
CALL  IMTPFT 
JO  PFTURN 
II  InO 

CALL  V-'rr|TAti»n1innti  rnis  we£fa«cm  $UPJ>D1?TC1}  I N  taut  m 

call  APTTATi?intft*ii]*«nT  tht  ur-jviRvjTr  fiF  LJ7AH  riHrLjTrp  *i 

CALL  ■  «  I  T  A  -  ,  ,  A  n  A  ,  *  sr  r  r  VC  IT  DimniF  Attn  THE  AnVirren  ‘I 

CALL  fTMf  flPiniShftO.tPF^'AHCK  MH1JFC  rr  A  GE  N-C  t  cr  T H fr  *1 
CALL  "PMj  f  UF  f 4  ?  A  NO  »t 

CALL  MM  T  ?  ( JTfl  •  *1  n Pj J  inrico  #V  rn*.  onA-C  MR  OEvriDPIff  rrT  n 

CALL  f  p  I  70  ,  ftnf? ,  I  CflTEP  ,  CpAFn  ,  Itrv  vriRr  |3#H0f  uUPE"  *J 

CALL  Alrt  TA74  ,  I  cn*i  TPArT  AT  )  tH  ft  Pi  | ?7  7  ,  f  | 

CAM.  r  Fjj 
C.i.L  i!>TP|M 
F  MO 

(31  FOR  PIC 

SlinROHTJMC  PICIF  ,K  ,f)F,  I  ) 

C  THIS  DISPLAYS  DIAGRAM  USED  IN  MAKING  MOVIE# 
i riTEGf.p  r  ci  non  > 
cai.i.  p  c  t  ,  2.  r>  a  ♦  r»  s  n  j 
CALL  Pns  (  2Sf.  ,S50  ) 

CALL  D  A  S  H ( 7  A  ft  t  s  S  0 ) 

CALL  PCS  (  ?Sft  ,Sf,0  ) 

CALL  DASM( 7Aft  ,9?‘M 
CALI.  POSIftl  *1.550) 

C  A  1. 1.  V  (  C  (  A  I  ft  ,  A  I  ?  ) 

CALL  POS(7ftn»SSO) 

CALL  V F  C  C  7  A  ft  #  9  ?  ft ) 

CALL  LfMAR 
00  l  Mr  I  0,370,1 2 
I  CALL  PnTipRA+MjSIS) 

CALI.  i¥RfTAT(J?l7t*J79,»  (I7HFYL  TO  DISPLAY  ■  ,  FM.I,|ftM  SCREEN  0  1  AM 
XRS  ) • , ( ) 

00  ?  Mr:  I  0,52"  ,  |  2 
?  CALL  OPT  l  >«*.♦••  ,3«A  ) 

CAI.I.  l"PlTAT(3tft,3SO,  •  MSHEYF.  TP  STAGE  »  ,IH,6M  FFET  l*,K) 

CALL  SFMPF 

nr turn 
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